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Abstract 


Experiments  have  been  conducted  in  the  GALCIT  Supersonic  Shear  Layer 
Facility  to  investigate  some  aspects  of  high-Reynolds-number,  turbulent,  shear- 
layer  flows  in  both  incompressible-  and  compressible-flow  regimes.  Experiments 
designed  to  address  several  issues  were  performed;  effects  of  inflow  boundary  con¬ 
ditions,  freestream  conditions  (supersonic/subsonic  flow),  and  compressibility,  on 
both  large-scale  dynamics  and  small-scale  mixing,  are  described. 

Chemically-reacting  and  non-reacting  flows  were  investigated,  the  former  rely¬ 
ing  on  the  (H2  +  NO)/F2  chemical  system,  in  the  fast-kinetic  regime,  to  infer  the 
structure  and  amount  of  molecular-scale  mixing  through  use  of  “flip”  experiments. 
A  variety  of  experimental  techniques,  including  a  color-schlieren  visualization  sys¬ 
tem  developed  as  part  of  this  work,  were  used  to  study  the  flows. 

Both  inflow  conditions  and  compressibility  axe  found  to  have  significant  effects 
on  the  flow.  In  particular,  inflow  conditions  are  “remembered”  for  long  distances 
downstream,  a  sensitivity  similar  to  that  observed  in  low-dimensionality,  non-linear 
(chaotic)  systems.  The  global  flowfields  (freestreams  coupled  by  the  shear  layer) 
of  transonic  flows  exhibit  a  sensitivity  to  imposed  boundary  conditions,  i.e.,  lo¬ 
cal  area  ratios.  A  previously-proposed  mode- selection  rule  for  turbulent-structure 
convection  speeds,  based  on  the  presence  of  a  lab-frame  subsonic  freestream,  was  ex¬ 
perimentally  demonstrated  to  be  incorrect.  Compressibility,  when  decoupled  from 
all  other  parameters,  c.g. ,  Reynolds  number,  velocity  and  density  ratios,  etc.,  re¬ 
duces  large-scale  entrainment  and  turbulent  growth,  but  slightly  enhances  small- 
scale  mixing,  with  an  associated  change  in  the  structure  of  the  molecularly-mixed 
fluid.  This  reduction  in  shear-layer  growth  rate  is  examined  and  a  new  parameter 
that  interprets  compressibility  as  an  energy-exchange  mechanism  is  proposed.  The 
parameter  reconciles  and  collapses  experimentally-observed  growth  rates. 
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CHAPTER  1 


Introduction  and  background 


During  the  past  three  decades,  turbulent  shear-layer  flows  have  been  the  subject 
of  rather  intense  experimental  and  numerical  scrutiny.  This  interest  is  due,  in  large 
part,  to  their  inherent  geometric  simplicity;  the  nearly-parallel,  isobaric  freestreams 
provide  a  spatially-uniform  boundary  condition  for  the  flow  (Fig.  1.1).  Despite 
this  simplicity,  the  response  of  the  flow  to  various  parameters  remains  a  matter  of 
ongoing  study  and  controversy.  On  a  practical  note,  as  a  relatively  efficient  mixer  of 
initially-separated  constituents,  this  flow,  and  others  like  it,  play  a  critical  role  in  the 
performance  of  many  industrial  devices,  ranging  from  chemical  lasers  to  envisaged 
hypersonic-propulsion  engines. 


Fig.  1.1  Color-schlieren  photograph  of  a  subsonic,  turbulent,  shear  layer  flow.  The 
flow  in  this,  and  all  subsequent,  flow  visualizations  is  from  left  to  right. 
The  subscript  ()i  refers  to  the  properties  of  the  high-speed  freestream  flow, 
and  the  subscript  ()2  refers  to  those  of  the  low-speed  freestream  flow. 
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This  thesis  explores  several  aspects  of  high-Reynolds-number,  turbulent,  shear- 
layer  flows,  extending  previous  work  in  both  the  incompressible-  and  compressible- 
flow  regimes.  In  particular,  experiments  to  investigate  effects  of  inflow  conditions, 
freestream  boundary  conditions,  and  compressibility  will  be  presented,  along  with 
flow-modeling  considerations  describing  these,  and  previous,  results.  Comprehen¬ 
sive  review  articles  discussing  turbulent  shear-layer  mixing  and  combustion  (Di- 
motakis  1991a),  incompressible-flow  combustion  (Coats  1996),  and  compressibility 
effects  on  turbulence  (Lele  1994)  are  relevant.  A  brief  overview  of  the  particular 
subjects  of  interest  is  offered  below. 


An  important  catalyst  for  this  recent  activity  was  the  discovery  of  large-scale, 
spanwise-coherent  structures  in  fully-turbulent,  incompressible,  shear-layer  flows 
(Brown  &  Roshko  1974).  Examples  of  these  structures  are  evident  in  Fig.  1.1. 
In  addition  to  casting  high-Reynolds-number  free-shear  flows  in  a  new  light,  this 
finding  spawned  a  series  of  models  that  recognize  the  importance  of  the  large-scale- 
structure  (Galilean)  rest  frame,  e.g .,  Brown  (1974)  and  Dimotakis  (1986).  These 
models,  extending  previous  work  for  uniform-density  flow,  e.g.,  Abramovich  (1963) 
and  Sabin  (1965),  exhibit  a  dependence  of  the  incompressible-flow  (M  =  0)*  growth 
rate,  8q(x)/x,  on  freestream-velocity  and  -density  ratios,  i. e., 


So 

x 


So 

X 


_  U2  _  P?.\ 

~U^S~  pj- 


(1.1) 


In  particular,  the  geometric  and  similarity  considerations  employed  by  Dimotakis 
(1986)  yield  an  estimate  for  a  spatially-growing  layer  given  by, 


— (r,s)  ~  C6 

X 


(1  —  r)(l  +  s1/2)  J,  (l-^Vfl  +  s1/2) 


2(1  +  51/2r) 
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+  2.9  (1  +  r)/(l 


l3_\ 

~r)j  ’ 


(1.2) 


where  C$  is  independent  of  r  and  s.  The  temporal-growth  model  of  Brown  (1974) 
provides  an  expression  equivalent  to  Eq.  1.2  above,  with  the  exception  of  the  term  in 
braces  that  accounts  for  the  upstream/downstream  asymmetries  arising  from  spatial 
growth.  Although  experimental  data  are  in  accord  with  the  overall  dependence 
suggested  by  Eq.  1.2,  there  is  substantial  variance  in  inferred  values  for  C$,  i.e., 


0-25  <  C6  <  0.45, 


(1.3) 


* 


The  subscript  ()o  denotes  incompressible  flow,  e.g Sq/x  is  the  shear-layer  growth  rate  at 
incompressible-flow  conditions. 
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especially  between  experiments  performed  in  different  facilities  (Dimotakis  1991a). 
Significantly,  initial,  or  inflow,  conditions  are  known  to  affect  far-field  properties  of 
the  turbulence,  including  growth  rates.  As  a  consequence,  it  is  widely  speculated 
that  such  an  influence  is  at  least  partially  responsible  for  the  range  exhibited  in 
Eq.  1.3.  Inflow  conditions  and  their  influence  on  shear-layer  growth  and  mixing  are 
discussed  in  Ch.  3. 

The  experiments  of  Brown  &  Roshko  (1974)  also  addressed  compressibility  ef¬ 
fects,  albeit  in  an  indirect  manner.  At  the  time  of  those  experiments,  it  was  well 
appreciated  that  Mach  number  had  a  significant  effect  on  shear-layer  growth.  In 
particular,  increasing  Mach  number  was  associated  with  decreasing  growth  rate. 
A  proposed  explanation  was  that  an  appropriate  transformation,  e.g.,  Howaxth- 
Doridnitsyn  type,  would  account  for  this  decreased  growth.  Such  an  explanation 
implicitly  assumes  that  compressibility  acts  by  modifying  the  density  field  in  a 
passive  manner.  Similar  explanations  have  been  relatively  successful  in  describing 
wall-bounded  flows,  up  to  Mach  numbers  of  5  or  so  (e.g.,  Bradshaw  1977).  The 
variable-density  experiments  of  Brown  &  Roshko,  however,  showed  that  the  incom¬ 
pressible  shear-layer  growth-rate  decreases  by  about  30%  for  $  —  1/7  (relative  to 
s  =  1),  much  less  than  the  300%  decrease  observed  for  an  equivalent  (M  ~  5.7) 
compressible  flow.  This  result  demonstrates  that  compressibility,  beyond  its  ef¬ 
fects  on  the  density  field,  fundamentally  alters  large-scale  dynamics  in  free-shear 
turbulence. 


The  success  of  convective-frame  views  in  describing  incompressible  flows  in¬ 
spired  an  extension  of  those  ideas  to  the  compressible-flow  regime.  To  scale  com¬ 
pressibility  effects,  Papamoschou  &  Roshko  (1988)  suggested  the  convective  Mach 


numbers, 


Mci 


Ui-Uc 

oi 


Mc  2 


Uc  U2 

0.2 


(1.4) 


where  U{  are  the  freestream  speeds,  ai  the  corresponding  freestream  speeds  of  sound, 
and  Uc  the  large-scale  structure  convection  speed,  extending  ideas  employed  by 
Mack  (1969,  1975,  and  1984)  in  supersonic,  boundary-layer  transition.  A  single 
parameter  was  subsequently  suggested  by  Papamoschou  (1989),  termed  here  the 
(total)  convective  Mach  number, 


Mc  = 


AU 
0\  +  02 


(1.5) 
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where  AU  =  U\  —  ?72  is  the  freestream-speed  difference.  For  equal  freestream 
specific-heat  ratios  (71  =  72),  extension  of  the  Dimotakis  (1986)  stagnation-point 
recovery  model  to  isentropic,  compressible  flow  (Papamoschou  &  Roshko  1988)  pre¬ 
dicts  symmetric  convection,  i.e.,  Mcl  =  Mc2  =  Mc,  with  all  three  values  relatively 
close  for  the  possible  range  of  7;’s. 

There  are  some  complicating  difficulties  with  such  a  view  in  the  compressible- 
flow  regime.  Linear-stability  analyses  ( e.g .,  Gropengiesser  1970,  Ragab  &  Wu  1988, 
Sandham  &  Reynolds  1990,  and  Zhuang  et  al.  1990)  indicate  the  potential  for 
several  distinct  “modes”  of  instability  in  unbounded  flow,  namely  fast,  central,  and 
slow  modes.  This  behavior  differs  from  incompressible  flow,  where  a  single  (central) 
mode  is  dominant  and  has  been  hypothesized  to  be  a  linear  “embryo”  of  the  two- 
dimensional  coherent  structures  discussed  above  (Roshko  1976). 

Schlieren/shadowgraph  visualizations  of  supersonic,  axisymmetric,  turbulent 
jets  indicate  fast-mode  existence  (Lowson  &  Ollerhead  1968  and  Tam  1971),  as  well 
as  slow-  and  fast-mode  coexistence  (Oertel  1979),  in  the  upstream,  annular-shear- 
layer  region.  Lowson  &  Ollerhead  (1968)  remarked  that  the  angle  of  the  Mach 
waves  generated  by  their  axisymmetric  jet  indicated  a  large-scale  structure  convec¬ 
tion  speed,  Uc,  of  about  70-80%  of  the  jet-exit  velocity;  a  value  significantly  higher 
than  that  inferred  for  incompressible  jets.  Following  these  observations,  evidence 
of  slow  and  fast  modes  in  two-dimensional  shear-layer  flows  were  reported  by  Pa¬ 
pamoschou  (1989)  and  subsequently  by  Hall  et  al.  (1993)  and  Papamoschou  &  Bun- 
yajitradulya  (1996).  These  (two-dimensional)  shear-layer  data  appear  to  indicate, 
however,  existence  of  either  the  slow  or  the  fast  mode  in  a  given  flow.  Papamoschou 
(1989)  suggested  an  empirical  mode-selection  rule,  with  a  selection  criterion  based 
on  the  presence  of  subsonic  lab-frame  freestream  flow.  If  a  (lab-frame)  criterion 
for  mode  selection  is  valid,  or  alternatively,  if  multiple-mode  coexistence  occurs, 
the  fundamental  idea  of  a  purely  convective-frame  analysis  must  be  questioned.  A 
discussion  of  these  issues  is  contained  in  Ch.  4. 

Partially  because  of  these  difficulties,  reports  of  shear-layer  data  have  typi¬ 
cally  concentrated  on  behavior  as  a  function  of  the  total  convective  Mach  num¬ 
ber,  Mc.  This  average  compressibility  measure  lacks  the  rigor  of  incompressible- 
flow  models,  since  turbulent  features  do  not  convect  with  the  (central/isentropic) 
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speed  implicit  in  Eq.  1.5,  at  least  beyond  low  levels  of  compressibility.  Addition¬ 
ally,  flow- visualization  data,  e.g .,  Clemens  &  Mungal  (1995),  are  indicative  of  a 
relative  reduction  in  well-organized,  two-dimensional,  coherent  structure  with  in¬ 
creasing  compressibility.  Nevertheless,  Mc,  as  well  as  the  “relative  Mach  number” 
(Mrei  =  2 Mc),  has  been  widely  used  in  assessing  compressibility  effects  on,  for 
example,  growth  rate.  In  particular,  compressibility  is  usually  assumed  to  act  in¬ 
dependently  of  freestream- velocity  and  -density  ratios,  yielding, 

-  a  j-(Mc)  ±  fn(r,s).  (1.6) 

Oo/X  0  0 

Support  for  this  type  of  scaling  can  also  be  found  in  the  linear- stability  analyses 
discussed  above.  There  have  been,  however,  experimentally-observed  violations  of 
Eq.  1.6,  e.g .,  in  Hall  et  al.  (1993),  which  are  discussed  in  Ch.  6. 

The  growth  rate  measures  the  evolution  of  the  largest  (outer)  length  scale 
and  can  be  regarded  as  the  first  step  in  a  sequence  of  shear-layer  mixing  pro¬ 
cesses  (Dimotakis  1991a).  When  a  chemical  reaction  between  two  freestream  con¬ 
stituents  is  of  interest,  e.g.,  as  in  non-premixed  combustion,  then  mixing  processes 
at  molecular  scales  must  be  addressed,  since  the  (molecular-interaction)  chemical 
processes/reactions  occur  at  those  scales.  In  particular,  molecularly-mixed  fluid  is 
the  end  result  of  a  multi-stage  process  that  can  be  regarded  to  occur  in  a  sequence 
of  Lagrangian  steps,  e.g.,  “stirring”  then  “mixing”  (Danckwerts  1958).  In  such  a 
process,  parameters  known  to  influence  the  largest  scales  of  the  flow,  i.e.,  inflow 
conditions,  velocity  and  density  ratios,  and  Mach  number,  must  be  augmented  by 
variables  that  influence  smaller  spatial  scales  in  the  flow,  e.g.,  Schmidt  number 
and  (local)  Reynolds  number.  Shear-layer  mixing  data  indicate  that  even  past  the 
mixing  transition,  where  a  qualitative  change  in  (small-scale)  mixing  occurs,**  i.e., 

Ree  s  >  104,  (1.7) 

there  is  a  dependence  of  molecular  mixing  on  Reynolds  number.  In  particular, 
both  experiments  (Mungal  et  ah  1985,  Frieler  &;  Dimotakis  1988,  and  Bond  & 
Dimotakis  1996),  as  well  as  models  (Broadwell  &  Briedenthal  1982,  Broadwell  & 


** 


A  brief  discussion  of  this  transition  can  be  found  in  Ch.  3.  For  a  more  complete  discussion,  see 
Dimotakis  (1993). 
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Mungal  1991,  and  Dimotakis  1987)  indicate  a  modest  decrease  in  molecularly-mixed 
fluid,  within  the  incompressible  shear  layer,  as  Reynolds  number  increases.  Such  a 
dependence  suggests  that  if  compressibility  effects  on  molecular  mixing  are  to  be 
ascertained  with  any  confidence,  variations  in  Reynolds  number,  as  well  as  other 
flow  parameters,  must  be  minimized.  Experiments  that  observe  this  restriction  are 
described  in  Ch.  5. 
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CHAPTER  2 

Facility,  instrumentation,  and  measurements 


2.1  Facility  overview 

The  GALCIT  Supersonic  Shear  Layer  (S3L)  Facility  is  a  two-stream,  blow¬ 
down,  wind  tunnel,  capable  of  attaining  freestream  Mach  numbers  spanning  from 
incompressible  to  highly-compressible  gas-phase  flow,  in  both  chemically-reacting 
and  non-reacting  configurations.  A  brief  overview  of  this  facility  is  presented  here 
for  completeness.  Additional  details  can  be  found  in  Hall  &  Dimotakis  (1989), 
Hall  (1991),  and  Bond  (1998).  Two  major  upgrades  undertaken  as  part  of  this 
work,  i.e.,  extension  to  bi-supersonic  flow  and  replacement  of  the  surge  tank,  will 
be  documented  here. 

The  facility,  shown  schematically  in  Fig.  2.1,  operates  with  a  run  time  of  trun  > 
tdat  —  2  —  4  s,  where  fdat  is  the  data-recording  time.  Gas  for  the  two  freestreams 
is  supplied  by  independent  (high-  and  low-pressure,  usually  corresponding  to  high- 
and  low-speed,  respectively)  flow  systems  charged  by  standard  bottled-gas  supplies 
prior  to  the  experimental  run.  This  facility  can  produce  flows  with  fast-kinetic 
reactants,  whose  chemical  reaction  and  associated  temperature  rise  yields  quanti¬ 
tative  measurements  of  molecular-scale  mixing  in  the  shear-layer  region.  In  such 
experiments,  a  dilute  mixture  of  H2  and  NO  is  carried  in  the  high-speed  stream 
and  a  dilute  concentration  of  F2  is  carried  in  the  low-speed  stream.  The  remaining 
diluent  components  are  selected  to  control  freestream  thermodynamic  properties, 
e.g. .  densities,  speeds-of-sound. 

To  provide  the  molecular  number  densities  required  for  “fast”  chemical  reac¬ 
tions,  large  (unit)  Reynolds  numbers,  and  high  optical-image  cross  sections,  the 
test  section  operates  at  nominally-atmospheric  static  pressure.  This  is  different 
from  most  other  supersonic-flow  facilities,  in  which  test  sections  operate  at  lower- 
pressure  conditions  to  ease  demands  on  upstream  flow-supply  systems. 
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Legend 

High-pressure  gas  flow  (H2+NO  leg) 
^  Low-pressure  gas  flow  (F2  leg) 
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Fig.  2.1  Overall  facility  gas-flow  schematic. 

For  the  shear-layer  flows  investigated  in  these  experiments,  the  high-pressure 
freestream  speed  is  maintained  by  plenum-pressure  control,  operated  in  either  an 
open-loop  (program-control)  configuration  at  low  mass-flux  conditions,  e.g.,  flows 
discussed  in  Ch.  3,  or  a  closed-loop  (nonlinear-feedback-control)  configuration  at 
moderate-to-high  mass-flux  conditions,  e.g.,  flows  discussed  in  Ch.  4.  Gas  is  supplied 
to  this  plenum  from  a  1.2  m3  pressure  vessel,  at  a  maximum  pressure  of  104kPa. 
During  blowdown,  vessel  pressure  can  drop  by  as  much  as  50%,  necessitating  an 
active-control  system.  The  gas  in  the  vessel  is  maintained  at  a  nearly  constant 
temperature,  however,  by  rolled  aluminum  screen  (large  heat  capacity)  that  occupies 
the  pressure-vessel  volume. 

The  mass  flux  through  the  low-pressure  plenum  is  set  by  a  choked  metering 
valve  kept  fixed  during  each  run.  Gas  to  this  plenum  is  supplied  from  a  0.57  m3 
vessel  that  encloses  a  bladder  bag  containing  the  F2-laden  charge.  This  vessel  is 
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pressurized,  in  turn,  by  a  N2-filled,  large-volume  surge  tank  whose  pressure  falls  by 
3  —  5%  during  a  run.  This  nearly-isobaric  arrangement  obviates  the  need  for  low- 
pressure-side  active  control,  greatly  simplifying  the  flow-control  process,  as  well  as 
providing  a  relatively  safe  (teflon-buffered)  method  for  handling  extremely-corrosive 
fluorine  gas. 


FIG.  2.2  Northeast  view  of  the  upgraded  surge  tank  installation. 


One  major  engineering  upgrade  undertaken  during  this  work  replaced  the  orig¬ 
inal  (circa  1940)  12.7  m3  surge  tank  and  associated  delivery  piping  with  a  higher- 
capacity  system.  The  original  system  had  a  maximum  working  pressure  of  700  kPa, 
limited  by  the  tank  pressure  rating.  This  significant  modification  was  necessitated 
by  the  planned  demolition  of  the  GALCIT  10-Foot  Wind  Tunnel  during  the  sum¬ 
mer  of  1997,  since  the  original  surge  tank  resided  adjacent  to  this  tunnel,  at  both 
basement  and  sub-basement  levels,  in  the  center  of  the  Guggenheim  Building.  The 
replacement  surge  tank,  a  custom  14.2  m3  unit  with  a  (maximum)  working  pressure 
of  1150  kPa,  was  installed  in  the  Karman/ Guggenheim  courtyard  (Fig.  2.2)  and  was 
connected  to  the  low-pressure  supply  system  in  July  of  1997.  The  upgraded  system 
has  since  been  successfully  tested  over  a  range  of  pressures  and  flowrates. 


Fig.  2.3  Photograph  of  shear-layer  flow-management  and  test-section  regions.  Gas 
flows  from  left  to  right.  Image  is  approximately  2m  in  length. 

A  photograph  of  the  flow-management  and  test-section  regions  is  shown  in 
Fig.  2.3.  Coarse-mesh-screen,  honeycomb,  and  a  series  of  fine-mesh-screen  sections 
are  used  for  turbulence  management  in  the  region  upstream  of  the  contractions. 
The  two  streams  are  separated  in  their  respective  plena  by  a  splitter  plate,  which  is 
flat  on  the  top  (high-speed)  surface  and  curved  on  the  bottom  (low-speed)  surface, 
as  shown  in  Fig.  2.3.  The  splitter  plate  has  a  trailing-edge  thickness  of  0.13  mm,  an 
included  angle  of  5°,  and  is  estimated  (using  General  Electric  gauges)  to  have  an 
average  surface  roughness  of  0.4  yum. 

For  subsonic  freestream  flows,  the  modular  “nozzle  blocks”  are  configured  as 
specially-machined  contraction  sections.  Their  shapes  are  described  by  6th -degree 
polynomials,  which  are  iteratively  optimized  using  a  potential-flow  solver  (Pepin 
&  Dimotakis  1989)  to  obtain  the  inviscid  outer-flow  solution,  and  the  method  of 
Thwaites  (e.g.,  White  1974)  to  obtain  the  boundary-layer  thickness(es).  The  op¬ 
timization  endeavors  to  maximize  flow  uniformity  and  to  minimize  boundary-layer 
momentum  thicknesses  at  the  contraction  exit,  while  ensuring  suppression  of  any 
flow  separation  and  three-dimensional  Taylor- Gortler  instabilities.  This  flow  in¬ 
stability  can  occur  in  the  boundary  layers  on  the  concave  surfaces  of  both  nozzle 
walls,  as  well  as  on  the  concave  part  of  the  lower  splitter-plate  surface.  The  design- 
iteration  process  minimizes  the  potential  for  this  instability  in  all  susceptible  regions. 
For  supersonic  freestream  flows,  these  nozzle  blocks  are  configured  as  converging- 
diverging  (Laval)  nozzles,  which  employ  similar  contraction  sections  to  the  throat, 
mated  to  supersonic-expansion  contours  designed  using  the  characteristic-capturing 
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method  of  Sivells  (1978).  As  seen  in  Fig.  2.3,  the  5°  included  angle  of  the  splitter 
plate  is  matched  on  the  bottom  nozzle  block,  i.e.,  the  bottom-freestream  flow  en¬ 
ters  the  test  section  at  a  5°  angle  relative  to  the  top-freestream  flow.  Flow  in  the 
test  section,  downstream  of  the  contractions/nozzles,  is  confined  by  top  and  bottom 
guidewalls,  with  optical  windows  serving  as  sidewalls,  with  a  span,  b  =  15.3  cm.  In 
these  experiments,  the  upper  guidewall  is  maintained  horizontal  (parallel  to  the  top 
splitter-plate  surface),  with  the  lower  guidewall  adjusted,  as  necessary  in  each  flow, 
to  remove  any  mean  streamwise  pressure  gradients. 

Upon  exiting  the  test  section,  the  flow  encounters  a  shower-tunnel  and  catch- 
bag  containment,  handling,  and  neutralization  system.  A  set  of  shower  heads  inside 
the  tunnel  spray  a  H20/NaOH  mixture  into  the  exiting  flow,  simultaneously  neu¬ 
tralizing  and  cooling  the  exhaust  gases.  The  neutralization  process  transforms  the 
exhaust  into  a  mixture  of  non-toxic  gases,  more  H20,  and  NaF,  the  active  ingre¬ 
dient  in  many  toothpastes.  The  150  m3 -volume,  atmospheric  pressure  catch  bag 
collects  and  contains  the  exhaust  mixture  downstream  of  the  shower  tunnel.  An 
additional  set  of  shower  heads,  within  the  bag,  provides  further  neutralization  and 
cooling  of  the  exhaust  products.  Finally,  in  preparation  for  eventual  venting,  the 
catch  bag  allows  for  controlled  dilution  of  any  potentially-flammable  gas  mixtures 
to  concentrations  comfortably  below  ignition  limits. 


2.2  The  upgrade  to  bi-supersonic  flow 

The  original  low-pressure  freestream  of  the  S3L  facility  was  limited  in  Mach 
number  to  subsonic  flow  (M2  <  0.5)  by  the  existing  components.  Consequently, 
an  upgrade  of  these  components  was  required  to  explore  the  regime  of  bi-supersonic 
flow,  where  both  freestreams  are  supersonic  in  the  laboratory  frame.  Two  major 
modifications  were  required  to  allow  bi-supersonic  operation  of  the  facility:  a  new 
control/regulation  system,  and  a  new  nozzle  block  designed  for  an  appropriate  Mach 
number.  Fortunately,  the  existing  low-pressure  side  reservoirs  (supply  and  surge 
tanks)  were  deemed  capable  of  supplying  enough  gas  for  adequate  run  times  at 
supersonic  Mach  numbers. 
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2.2.1  Low-pressure  freestream  control  system 

One  might  expect  that  an  active-control  valve,  similar  to  the  unit  which  reg¬ 
ulates  the  high-pressure  freestream  flow,  would  be  required  to  compensate  for  the 
anticipated  stagnation-pressure  drop  during  a  run  (35  kPa  from  the  original  surge- 
tank  pressure  of  700 kPa).  However,  since  the  existing  high-pressure  control  valve 
can  be  pre-programmed,  it  is  possible  to  achieve  the  desired  pressure  ratio  across 
the  splitter-plate  tip  throughout  a  run  with  fixed-area- valve  metering  of  the  low- 
pressure  freestream.  In  order  to  accommodate  increased  mass  fluxes  associated  with 
supersonic  Mach  numbers,  an  increase  in  (sonic)  valve-orifice  area  was  required.  To 
provide  this  increase,  a  new  sleeve  was  machined  by  the  Aeronautics  Machine  Shop 
with  approximately  triple  the  choked-flow  area  of  the  original  (subsonic)  sleeve. 
These  modular  sleeves  are  easily  interchanged,  increasing  the  operating  envelope  of 
the  S3L  facility  to  include  bi-supersonic  flow,  in  addition  to  high-mass-flux  subsonic- 
flow  conditions,  e.g.,  flows  discussed  in  Ch.  5. 

The  existing  passive-blowdown  method  of  metering  the  low-pressure  freestream 
has  several  significant  advantages  over  an  active-control  scheme.  In  particular, 
this  allows  continued  use  of  the  original  operating  procedure  and  control  algo¬ 
rithm/computer,  minimizing  changes  to  a  well-tested  facility.  Additionally,  all 
existing  low-pressure  “settings”  remain  valid,  permitting  low  mass-flux  runs  to  oc¬ 
cur  with  a  minimum  of  recalibration.  Finally,  omission  of  a  complex  valve  assembly 
eliminates  the  significant  safety  considerations  of  long-term  fluorine  passivation  and 
compat  ability. 
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2.2.2  Supersonic  nozzle 

A  preliminary  study  indicated  that  a  Mach  number  of  less  than  approximately 
1.13  would  be  achievable  using  existing  reservoir  components.  This  limit  is  dictated 
if  a  minimum  allowable  run  time  of  2  s  is  desired,  e.g.,  to  provide  adequate  data- 
acquisition  record  lengths,  and  a  nozzle-exit  height  of  38mm  is  desired,  e.g.,  to 
admit  the  shear-layer  growth  expected  over  the  test-section  length.  Near-sonic  Mach 
numbers,  however,  dictate  throat-to-expansion  area  ratios  correspondingly  close  to 
unity.  Considering  uncertainties  inherent  in  the  displacement-thickness  estimate^  of 
0.35  mm,  as  well  as  errors  incurred  in  machining  a  large  block  of  aluminum,  it  seemed 
wise  to  design  a  nozzle  with  a  throat-to-exit  expansion  as  large  as  possible,  i.e.,  a 
nozzle  with  as  high  a  Mach  number  allowed  by  the  above  constraints.  Consequently, 
a  Mach  1.13  nozzle,  with  a  throat-to-exit  expansion  of  approximately  0.85  mm,  was 
fabricated  by  Caltech  Central  Engineering  Services  and  was  used  throughout  the 
bi-supersonic  flow  experiments  described  in  Ch.  4. 


2.3  Method,  instrumentation,  and  diagnostics 

Three  types  of  experimental  results  will  be  documented  here:  schlieren  visu¬ 
alizations  (both  black-and-white  and  color),  mean-temperature  measurements,  and 
mean-pressure  measurements.  Schlieren  visualization  records  instantaneous,  overall, 
spanwise-averaged  (two-dimensional)  structure  in  the  layer,  indicating  whether  the 
shear-layer  growth  is  linear  and  also,  depending  on  the  selected  sensitivity,  whether 
any  significant  acoustic  disturbances  are  present.  Additionally,  the  visual-thickness 
growth  rate,  6V 5s/x,  can  be  measured  from  such  data.  When  configured  correctly, 
color  schlieren  registers  the  additional  dimension  of  information  provided  by  color 
to  record  both  components  of  line-of-sight-averaged  index-of-refraction  gradients  in 
the  flow  (see  Settles  1985  and  Appendix  A),  as  opposed  to  the  single  component 
recorded  by  traditional,  black-and-white,  knife-edge  systems. 


t  Estimate  provided  by  Spalding-Chi  turbulent  boundary-layer  model  internal  to  the  Sivells 
nozzle- design  code. 
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The  temperature  profiles  are  recorded  by  a  rake  of  sixteen,  50  ^m-diameter, 
chromel/alumel,  exposed-junction  thermocouples,  that  produce  a  running  temporal 
average  at  the  respective  measurement  locations.  These  probes  have  a  response 
time  estimated  to  be  approximately  0.02  s.  An  adjacent  rake  of  sixteen,  1.5  mm- 
diameter  tubes  record  pitot  pressures,  with  static  pressures  along  both  top  and 
bottom  guidewalls  monitored  by  flush-mounted  transducers.  As  with  temperature 
data,  pressure  data  represent  running  temporal  averages  at  the  various  measurement 
stations.  Both  temperature  and  pressure  measurements  are  recorded  within  3  mm 
of  the  midspan  plane,  with  the  rakes  positioned  near  the  test-section  exit  (Fig.  2.3) 
to  minimize  probe-induced  flow  blockage. 

These  data  are  processed  by  a  set  of  analog  3-pole  Butterworth  amplifier/filter 
banks  with  knee  (—3  dB)  frequencies  set  at  100  Hz.  These  are  sampled  at  200  Hz 
by  a  National  Instruments  analog-to-digital  converter,  operating  under  control  of 
a  Gateway-PC/LabView  environment  developed  for  these  experiments.  The  data 
are  post-processed  by  truncating  the  flow  startup  and  shutdown  periods,  to  obtain 
the  steady-flow  period  of  the  run.  This  record  is  then  divided  into  eight  segments, 
each  of  which  is  individually  processed  to  check  for  discrepancies  with  respect  to  the 
entire  ensemble.  An  example  of  such  an  ensemble  of  (eighth-segment)  temperature 
profiles  is  presented  in  Fig.  2.4.  In  this  example,  the  maximum  standard  devia¬ 
tion  of  the  mean  (of  16  independent  thermocouple  records)  is  less  than  2%  of  the 
mixture  adiabatic-flame  temperature  rise,  ATf.  The  inferred  mixed-fluid  quantities 
described  below,  e.g.,  Eq.  2.6,  are  based  on  the  complete  steady-flow  ensemble,  and 
are  subject  to  a  statistical  uncertainty  closer  to  1%  of  ATf. 

Non-reacting  flow  temperatures  are  used  to  construct  a  “cold” -temperature 
profile  for  each  flow  condition.  This  profile  serves  as  a  reference,  and  when  sub¬ 
tracted  from  reacting-flow  mean-temperature  measurements,  yields  the  temperature 
rise,  AT,-  =  A T(x,y,),  at  each  measurement  location.  The  resulting  temperature 
rises  are  divided  by  the  adiabatic  flame  temperature  rise  of  the  reacting  mixture, 
ATf,  to  yield  normalized  mean  temperature-rise  measurements,  AT, /ATf.  These 
are  least-squares  fit  with  a  functional  form, 

@(y)  =  =  exP  [-  (po  +  pi  y  +  P2  y2  +  pz y3  +  p 4  y4 )]  ,  (2.1) 
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Fig.  2.4  Eight  records  of  raw  temperature  data.  Flow  is  (tripped)  Case  3.1  condi¬ 
tions  detailed  in  Ch.  3,  below. 

which  is  then  used  to  estimate  the  quantities  described  below. 

For  incompressible-flow  experiments,  e.y.,  flows  discussed  in  Ch.  3,  velocity  data 
were  deduced  from  direct  pitot-  and  static-pressure  measurements,  as  described  in 
Hermanson  &  Dimotakis  (1989), 

u(v)  -  (2'2) 

where  Pp(y )  is  the  pitot-pressure  measurement,  Pg  is  the  gui dewall  (static)  pressure 
measurement  and  p(y)  is  the  density.  Local  values  of  density  axe  estimated  as, 

p(y)  _  To 

Po  T0  +  AT(y), 


(2.3) 
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using  the  average  temperature  rise  measured  by  the  thermocouples,  as  well  as  the 
ambient  temperature  and  density,  To  and  po ,  respectively.  Velocity  profiles  are  then 
least-squares  fit  with  an  analytical  form, 

U{y)  =  A  +  B  tanh  (g0  +  qi  y  +  y2  +  93  y3  )  •  (2.4) 

For  compressible-flow  experiments,  pitot-pressure  measurements  were  only  used  to 
measure  high-  and  low-speed  freestream  velocities,  U\  and  U2.  respectively. 


As  alluded  to  above,  the  molecularly-mixed  fluid  was  sampled  using  kinetically- 
fast  chemical  reactions  to  perform  a  “flip”  experiment,  i.e.,  chemically- reacting  flows 
at  high  Damkohler  number, 

Da  =  >  1  ,  (2.5) 

7"chem 

where  rmjx  is  the  fluid-mechanical,  molecular-mixing  time  and  rchem  is  the  chemical- 
reaction  (completion)  time.  A  detailed  summary  of  the  rationale  and  measures 
employed  in  such  an  experiment  is  found  in  Dimotakis  (1991a)  and  Bond  (1998). 
A  brief  discussion  is  offered  below. 

In  this  high-Da  regime,  chemical-product  formation  and  heat  release  (temper¬ 
ature  rise)  provide  a  resolution-insensitive  probe  of  molecular  mixing,  mitigating 
the  spatial-sampling/resolution  difficulties  of  direct-measurement  studies  (c/.  dis¬ 
cussion  in  Breidenthal  1981,  Koochesfahani  &  Dimotakis  1986,  Dimotakis  1991a, 
and  Karasso  &  Mungal  1996).  In  order  to  attain  this  fast-chemistry  limit,  in  high¬ 
speed  flows,  use  of  the  (H2  +  N0)/F2  chemical  system  is  required. 


If  the  flow  is  of  gases  with  uniform  heat  capacity,  the  normalization  of  the  result¬ 
ing  temperature-rise  profile,  AT(y),  at  the  measuring  station,  x,  by  the  adiabatic 
flame  temperature  rise,  ATf,  as  in  Eq.  2.1,  yields  a  dimensionless  chemical-product 
profile.  As  chemical-product  formation  can  only  take  place  where  entrained  fluids 
are  mixed,  this  can  be  used  to  define  the  1%  shear-layer  (temperature-rise)  thick¬ 
ness,  <5t,  as  the  span  between  the  two  points  in  the  profile  at  1%  of  the  peak  in 
0(y).  This  has  been  found  to  be  in  good  agreement  with  alternative  measures  of 
the  local  shear-layer  width,  such  as  the  visual  thickness,  £vjs.  This  profile  is  also 
used  to  compute  the  chemical-product  thickness, 


/OO 

@(y) dy  > 

-00 


(2-6) 
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which  represents  the  thickness  of  an  equivalent  layer  that  has  uniformly-risen  to  the 
adiabatic-flame  temperature.  Finally,  the  product  thickness  is  normalized  by  the  1% 
temperature-rise  profile  thickness,  <$t,  to  obtain  the  dimensionless  product  fraction, 
(5^  J <5t  ,  providing  a  measure  of  chemical- product  formation  that  is  insensitive  to 
changes  in  the  local  shear-layer  width. 


Given  product  fractions  at  stoichiometric-mixture  ratios,  <j>  =  <f>o  and  <f>  =  1/ 
for  large  <j>0,  the  mixed-fluid  fraction  in  the  fast-kinetic  regime  can  be  estimated 
from, 


1  €(f,Q 


-  J  J  **.*)«*»  “I+I7S 


-oo  e*0 


(2.7a) 


where, 


7  = 


JL 


e4>o  - 


e.g. 


6s  a  0.056 


(2.7b) 


2  (1  4-  <f>o) 

In  this  expression,  X  is  the  high-speed  freestream-fluid  mole  fraction,  and  p(X,rj) 
is  its  probability  density  function  (pdf)  at  the  normalized  transverse  coordinate, 
rj  =  y/6r . 


For  a  flow  in  the  kinetically-fast  regime,  where  chemical-product  formation  is 
mixing-rate-limited  rather  than  kinetic-rate- limited,  Sm/ serves  as  a  resolution- 
insensitive  measure  of  the  mole  fraction  of  molecularly-mixed  fluid.  Operationally, 
attainment  of  fast  chemistry  (high-Da)  can  be  ascertained  by  demonstrating  that 
the  amount  of  chemical  product  remains  unaltered  as  the  chemical-kinetic  rate  is 
varied.  This  test  was  performed  by  proportionally  changing  the  concentrations  of 
all  reactants,  while  holding  <f>0  constant.  This  changes  the  kinetic  rate  by  not  only 
changing  the  number  density  of  reactant  species,  but  also  the  temperature  and 
associated  Arrhenius  kinetics  constants,  resulting  in  an  increase/decrease  of  Tchem) 
while  holding  constant. 


Such  measurements  also  yield  the  relative  amounts  of  entrained  freestream 
fluids  that  are  molecularly  mixed.  In  particular,  the  mixed-fluid  composition  (mole) 
ratio  can  be  calculated  as, 


Sn  = 


(V*r) 


1/00 


(2.8) 
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In  the  limit  of  large  <f>, o,  this  approaches  the  number  of  moles  of  high-speed  freestream 
fluid,  per  mole  of  low-speed  freestream  fluid,  that  are  molecularly  mixed  within  the 
shear-layer  region. 

The  predominant  mixed-fluid  composition,  Xe,  can  be  expected  to  reflect  the 
mole  ratio,  En,  of  fluid  entrained  from  the  two  freestreams  at  large  scales  (Dimotakis 
1986,  1991b),  i.e., 

<2-9> 

The  mixed-fluid  composition  ratio,  £n,  as  estimated  above,  however,  will  reflect  the 
overall  mixed-fluid  composition  and  need  not  match  En,  as  also  indicated  by  various 
mixing  models,  e.g.,  Broadwell  &  Breidenthal  (1982),  Broadwell  &  Mungal  (1991), 
and  Dimotakis  (1987). 

Effects  of  heat  release  on  the  flow  were  addressed  by  Hermanson  &  Dimo¬ 
takis  (1989).  As  was  noted  in  those  investigations,  the  most  significant  heat-release 
effect,  at  the  values  of  ATf/To  in  these  experiments,  is  the  dilatation-induced  (dis¬ 
placement)  flow  acceleration  in  the  confined-flow  environment.  In  the  present  ex¬ 
periments,  such  effects  were  accommodated  by  adjusting  the  lower  guidewall  to 
maintain  a  zero  streamwise  pressure-gradient  (no  acceleration),  for  all  cases,  ap¬ 
proximating  an  unbounded  shear-layer  flow,  as  noted  above.  Except  for  high  heat- 
release  conditions,  i.e.,  ATf/To  >  2  or  so,  only  small  effects  were  observed  on  the 
mixed-fluid  quantities.  As  noted  by  those  authors,  as  well  as  Dimotakis  (1991a), 
the  adiabatic-flame-temperature  rise  is  achieved  in  a  relatively  small  fraction  of  the 
fluid  within  the  layer,  with  an  average  temperature  rise,  distributed  throughout  the 
mixed  fluid,  that  is  significantly  lower.  This  can  can  be  seen  directly  in  the  mean 
temperature-rise  data  presented  later,  e.g.,  Fig.  3.1. 
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CHAPTER  3 


Effects  of  inflow  conditions 


This  chapter  details  results  from  a  set  of  incompressible-flow  experiments,  at 
high  Reynolds  number  (Reg  =  pAU  S(x)/p  ~  2  x  105),  in  which  inflow  conditions 
of  shear-layer  formation  were  varied.  These  variations  were  effected  by  perturb¬ 
ing  each,  or  both,  of  the  boundary  layers  on  the  splitter  plate  separating  the  two 
freestream  flows,  upstream  of  shear-layer  formation.  Results  of  chemically-reacting 
flip  experiments  reveal  that  seemingly  small  changes  in  inflow  conditions  can  have 
a  significant  influence  on  not  only  large-scale  structure  and  shear-layer  growth  rate, 
as  had  been  documented  previously,  but  also  on  chemical-product  formation  and 
molecular  mixing,  far  downstream  of  the  inflow  region. + 


$  These  results  were  reported  in:  Slessor,  Bond,  and  Dimotakis  (1998),  Turbulent  shear-layer 
mixing  at  high  Reynolds  numbers:  effects  of  inflow  conditions,  GALCIT  Report  FM98-1,  from 
which  this  Chapter  largely  derives. 
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3.1  Introduction 

It  is  an  implicit  assumption  of  turbulence  and  dynamic-similarity  theory  that 
for  high-enough  Reynolds  numbers  and  in  regions  of  the  flow  sufficiently  far  from 
inflow  boundaries,  the  behavior  of  the  flow  can  be  characterized  by  local-flow  pa¬ 
rameters,  at  least  in  a  statistical  sense  ( e.g .,  Batchelor  1953).  In  free-shear  layers 
and  some  other  flows,  however,  it  is  known  that  the  initial/inflow  conditions  of  the 
flow  can  exert  a  significant  influence  on  the  far-field  properties,  even  at  Reynolds 
numbers  that  are  considered  high. 

A  comparison  of  two  of  the  earliest  turbulent  “half-jet”  experiments,  by  Liep- 
mann  &  Laufer  (1947)  and  Wygnanski  &  Fiedler  (1970),  for  example,  revealed  a 
significantly-la.rger  growth  rate  in  the  latter  experiment,  a  discrepancy  eventually 
attributed  to  differences  in  the  flow  upstream  of  the  separation  (shear-layer  forma¬ 
tion)  region.  In  the  former  experiment,  the  boundary  layer  was  laminar,  while  in  the 
latter  it  was  turbulent.  This  increase  in  half-jet  growth  associated  with  a  turbulent 
splitter-plate  boundary-layer  initial  condition  was  noted  by  Batt  (1975),  and  con¬ 
firmed  by  Oster  et  al.  (1976).  Interestingly,  for  non-zero  velocity  ratios,  i.e.,  in  shear 
layers  with  a  coflowing  low-speed  stream,  Oster  and  coworkers  found  that  a  turbu¬ 
lent  high-speed  boundary  layer  had  the  opposite  effect,  i.e.,  a  reduced  downstream 
growth  rate.  This  behavior  for  co-flowing  freestreams  was  confirmed  by  the  subse¬ 
quent  work  of  Browand  &  Latigo  (1979),  Weisbrot  ei  al.  (1982),  and  Mungal  et  al. 
(1985).  In  much  of  this  work,  it  was  found  that  accompanying  velocity-fluctuation 
(rms)  and  Reynolds-stress  profiles  were  subject  to  significant  changes,  in  response 
to  perturbed  inflow  conditions. 


Shear  layers  are  not  unique  in  exhibiting  such  behavior.  Other  examples  of 
turbulent  flows  that  depend  on  initial/inflow  conditions  have  also  been  observed. 
Wygnanski  et  al.  (1986)  reported  far-field  properties  of  turbulent  wakes  that  de¬ 
pended  on  the  shape  of  the  body  that  produced  the  wake,  e.g.,  airfoil,  plate,  etc. 
A  sensitivity  to  initial  conditions  was  also  reported  for  wall-bounded  flows  by  Coles 
(1965),  in  transitional  Couette-flow  experiments,  a  behavior  corroborated  by  recent 
numerical  simulations  (Yao  &;  Moulic  1995). 
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A  requirement  that  the  local  Reynolds  number  must  be  appropriately  large 
stems  from  the  need  to  have  Reynolds  stresses  dominate  viscous  stresses,  but  also 
such  that  sufficient  small-scale  three-dimensionality  can  have  developed  and  be  sus¬ 
tainable.  In  shear  layers,  a  relatively-well-defined  transition  to  a  state  of  improved 
mixing  was  documented  by  Konrad  (1976),  in  gas-phase  shear  layers,  and  subse¬ 
quently  confirmed  in  liquid-phase  shear  layers  (Breidenthal  1978,  Bernal  et  al.  1979, 
Koochesfahani  &  Dimotakis  1986,  and  subsequently  others).  This  transition  was 
documented  to  occur  at  the  same  minimum  local  Reynolds  number,  based  on  the 
local  shear-layer  thickness  and  velocity  difference,  for  both  gas-  and  liquid-phase 
layers,  indicating  that  this  requires,  at  a  minimum, 

=  P± >  104  ,  (3.1) 

V- 

or  so.  Similar  transitions,  at  comparable  Reynolds  numbers,  have  been  documented 
in  many  other  flows,  so  this  behavior  and  the  applicability  of  this  criterion  must  be 
regarded  as  not  limited  to  shear-layers  ( cf  Dimotakis  1993). 

This  minimum  Reynolds  number  criterion  must,  however,  be  regarded  as  only 
necessary;  good  mixing  and  quasi-asymptotic  behavior  is  documented  to  have  ad¬ 
ditional  requirements.  Bradshaw  (1966)  indicated  that  shear-layer  turbulence  re¬ 
quired  a  relatively  long  distance  downstream  to  attain  an  asymptotic,  quasi-self- 
similar  state.  In  particular,  he  proposed  that  such  a  state  required  a  minimum 
downstream  distance,  i.e., 

x/e  >  103,  (3.2) 

where  x  is  the  distance  downstream  of  separation,  or  splitter-plate  trailing  edge, 
and  8  is  the  momentum  thickness  of  the  separating  boundary  layer;  an  inviscid 
criterion  that  can  be  varied  independently  of  Reynolds  number. 

It  was  noted  by  Dimotakis  &  Brown  (1976)  that  such  a  criterion  could  be  recast 
in  terms  of  a  minimum  number  of  large-scale-structure  interactions  ( e.g .,  pairings, 
or  triplings,  etc.),  estimated  by  them  to  be  approximately  four,  as  measured  in 
terms  of  an  estimated  number  of  mergings, 

m{x)  =  log2  >  4  ,  (3.3) 
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that  assumes  that  pairing  interactions  are  dominant,  i.e.,  a  measure  of  the  down¬ 
stream  distance,  x ,  as  scaled  by  the  initial-disturbance  spacing,  £$.  They  argued 
that  this  minimum  distance  and  minimum  number  of  interactions  were  required  to 
amortize  the  effects  of  initial  (inflow)  conditions. 


In  subsequent  work  by  Huang  &  Ho  (1990),  a  similar  scaling  of  the  downstream 
distance  was  proposed, 


1  —  r  /  x  \  ~  1  —  r  (  x  \ 
1  +  r  ~  1  +  r  \3O0i  / 


(3.4a) 


where  r  =  U2/U1  is  the  shear-layer  freestream  velocity  ratio  and  6\  is  the  initial, 
high-speed  boundary-layer  momentum  thickness.  Huang  &  Ho  argued  that  this 
parameter  effectively  described  the  approach  to  self-similarity  for  a  variety  of  shear- 
layer  flows,  i.e.,  for  different  velocity  ratios,  r  (cf.  also  Ho  &  Huerre  1984).  They 
suggested  that  a  minimum  value  of  8  was  required  for  self-similarity  of  the  peak 
turbulence  level  in  gas-phase  flows,  i.e., 


-Phh  >  8  , 


(3.4b) 


corresponding  to  the  location  of  the  second  large-structure  pairing,  according  to 
their  estimate,  i.e.,  a  criterion  more  relaxed  than  the  one  suggested  by  Dimotakis 
&  Brown.  These  authors  also  noted  that  the  properties  of  this  self-similar  state 
appeared  to  be  influenced,  if  not  determined  by,  the  initial  conditions  of  turbulence 
formation,  regardless  of  the  distance  from  separation. 


We  note  that  the  velocity-ratio  scaling  factor  in  Eq.  3.4  is  proportional  to  the 
local,  normalized,  large-scale  spacing,  £(x)/x,  independently  of  the  freestream  den¬ 
sity  ratio  (cf.  Dimotakis  1986),  and,  therefore,  as  a  measure,  the  Huang  Sz  Ho 
parameter,  P,  is  proportional  to  £(x)/£o,  rather  than  x/£q. 

Much  of  the  previous  work  included  measurements  of  shear-layer  growth  rates, 
Reynolds  stresses,  etc.  Only  limited  data  are  available,  however,  regarding  the 
response  and  sensitivity  of  the  shear-layer  concentration  field  and  its  mixing  prop¬ 
erties  to  changes  in  inflow  conditions.  Given  the  large-scale  effects  observed  in  the 
work  cited  above,  one  might  also  expect  significant  influences  in  the  small-scale  be¬ 
havior  of  the  flow,  which  is  important  in  molecular/diffusion-scale  mixing  processes. 
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Indeed,  hints  of  such  behavior  can  be  found  in  other  flows.  By  way  of  example,  far- 
field  jet-fluid  concentration-fluctuation  statistics,  in  the  far  field  of  a  turbulent  jet, 
were  found  to  be  strongly  influenced  by  small  (when  scaled  by  the  jet-exit  velocity) 
velocity  fluctuations  in  the  reservoir  fluid  that  the  jet  discharged  into  (Miller  1991). 

Some  recent  work  on  scalar  mixing  in  a  shear-layer  flow  was  reported  by  Karasso 
&  Mungal  (1996).  These  authors  found  the  Huang  t  Ho  (1990)  scaling  parame¬ 
ter  useful  in  their  analysis  of  their  measured  estimates  of  the  concentration-field 
probability-density  function  (pdf).  Dubbing  P  the  “pairing  parameter”  (Eq.  3.4), 
they  proposed  a  minimum  value  for  a  converged  mixed-fluid  state  (scalar  pdf)  of 
P  ~  22,  corresponding  to  a  distance  between  the  average  locations  of  third  and 
fourth  vortex  interaction;  more  closely  in  accord  with  the  original  proposal  of  Di- 
motakis  &  Brown  (1976).  Despite  expressing  reservations  about  the  accuracy  of 
their  direct-measurement,  laser-induced  fluorescence  results,  on  which  their  pdf 
experimental-estimates  relied,  they  suggested  that  this  value  determined  the  tran¬ 
sition  between  a  “non-marching”  pdf,  where  the  most-probable  concentration  is 
nearly-independent  of  the  transverse  coordinate,  y,  and  a  “tilted”  pdf,  in  which  the 
composition  is  correlated  with  the  local  mean,  i.e.,  a  function  of  y.  In  particular, 
they  proposed  that  the  former  would  be  observed  for  P  <  22,  while  the  latter  for 
P  >  22,  i.e .,  that  self-similarity  can  be  expected  for, 

Pkm  >  22  .  (3.4c) 

Karasso  &  Mungal  also  studied  the  effects  of  tripping  the  initial  boundary 
layers,  stating  that  the  effect  of  tripping  was  to  decrease  the  (absolute)  distance, 
x,  required  to  reach  self-similarity,  in  accord  with  the  behavior  reported  by  Bell 
&  Mehta  (1990),  who  stated  that,  “The  distance  required  to  achieve  self-similarity 
was  distinctly  shorter  in  the  tripped  case,  ...”  Karasso  &  Mungal  further  suggested 
that  beyond  this  distance,  inflow  conditions  were  “forgotten”  by  the  turbulence, 
with  a  unique  far-field  configuration,  that  only  depends  on  the  (local)  freestream- 
velocity  and  -density  ratios.  Since  tripping  the  boundary  layer(s)  can  be  expected 
to  result  in  an  increase  in  the  initial  momentum  thickness(es),  it  would  appear  that 
this  observation  is  not  in  accord  with  a  minimum-threshold  criterion  based  on  the 
pairing  parameter  P,  or  any  other  measure  of  distance  downstream  scaled  by  initial 
boundary-layer  thickness  scales. 
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3.2  Experiments 

Experiments  were  performed  at  a  single  nominal  flow  condition,  with  freestream 
speeds, 

Uo 

U\  ~  lOOm/s  and  U2  —  40m/s  ,  i.e.,  r  =  —  ~  0.4  ,  (3.5a) 

U 1 

and  a  density  ratio  of, 

s  =  —  ~  1.0  .  (3.5b) 

Pi 

The  velocity  and  density  ratios  were  chosen  to  facilitate  a  comparison  with  previous 
work  on  (incompressible)  chemically-reacting  flows,  at  similar  freestream  conditions 
( e.g .,  Mungal  et  al.  1985,  Frieler  &  Dimotakis  1988,  Hermanson  &  Dimotakis  1989, 
Bond  &  Dimotakis  1996),  as  well  as  Island  (1997),  who  investigated  molecular  mix¬ 
ing  in  shear-layers  at  different  conditions,  however,  using  a  fluorescence- quenching 
(“cold-chemistry”)  technique. 

The  pressure  and  temperature  profiles  in  these  experiments  were  recorded  at  a 
single  measuring  station,  x  =  0.365  m,  where  x  is  measured  from  the  location  of  the 
splitter-plate  tip.  Profiles  were  measured  as  a  function  of  the  transverse  coordinate, 
y,  that  is  also  referenced  to  the  splitter-plate  tip.  The  resulting  shear  layers  have  a 
local  Reynolds  number  at  the  measuring  station  in  the  range, 

1.8  x  105  <  Res  =  <  2.4  x  105  ,  (3.6) 

/X 

where  A 17  =  U1—U2  is  the  freestream  velocity  difference,  and  p  and  p  are  the  mixed- 
fluid  density  and  viscosity,  respectively.  These  are  calculated  using  the  properties  of 
a  gas  mixture  with  a  high-speed  fluid  mole  fraction  given  by  Xe  (Eq.  2.9),  assuming 
conservation  of  (total)  enthalpy  to  estimate  the  mixed-fluid  density,  p,  and  the 
method  of  Wilke  (1950)  to  estimate  the  mixed-fluid  (shear)  viscosity,  p.  These 
parameters  were  kept  nominally  fixed  to  isolate  effects  of  variable  inflow  conditions, 
as  each  is  known  to  influence  shear-layer  growth,  entrainment,  and  mixing.  The 
range  in  Eq.  3.6  arises  primarily  from  variations  in  local  shear-layer  thickness,  6(x ), 
in  each  case,  as  will  be  documented  below. 


High-speed  side: 
Low-speed  side: 


0/mm 

0.11 

0.17 


Ree  =  pU9/fi  x/6  P  (Eq.  3.4) 

640  3.3  x  103  47 

400  2.1  x  103 

Table  3.1:  Untripped  (natural)  splitter-plate  boundary-layer  flows. 

Changes  in  inflow  conditions  were  effected  through  various  combinations  of 
trip  wires  on  the  splitter  plate.  For  freestream  conditions  in  these  experiments, 
estimated  boundary-layer  properties  at  the  splitter-plate  tip  (x  =  0)  are  listed  in 
Table  3.1.  As  described  in  Ch.  2,  three  sets  of  experimental  runs  were  performed  to 
create  a  “flip  experiment”.  These  relied  on  chemically-reacting  flows,  with  stoichio¬ 
metric-mixture  ratios  of  <f>  =  (f>o  and  <j>  =  l/<po ;  for  Oo  =  8,  in  addition  to  non¬ 
reacting  flows  (“cold”  runs).  Table  3.2  lists  the  flow  conditions  and  freestream  fluid 
compositions  for  these  runs. 

3.1  -  Fluorine-rich  (<j>  =  8,  ATf  ~  171 K): 

lh  ~  lOOm/s  :  8.00%  Ar,  15.07%  He,  0.93%  H2,  0.14%  NO,  75.86%  N2 

P2  ~  40 m/s  :  16.00%  He,  8.00%  F2,  76.00%  N2  . 

3.2  -  Hydrogen-rich  ( <j>  —  1/8,  ATf  ~  267  K): 

Pf  ~  lOOm/s  :  8.00%  Ar,  4.00%  He,  12.00%  H2,  0.75%  NO,  75.25%  N2 

P2  —  40 m/s  :  6.45%  Ar,  16.00%  He,  1.55%  F2,  76.00%  N2  . 

3.3  -  Hydrogen-rich,  reduced  concentration  ( <f>  =  1/8,  ATf  ~  171 K): 

Ui  ~  lOOm/s  :  8.00%  Ar,  8.00%  He,  8.00%  H2,  0.50%  NO,  75.50%  N2 

U2  —  40m/s  :  6.97%  Ar,  16.00%  He,  1.03%  F2,  76.00%  N2  . 

3.4  -  Non-reacting: 

Pf  ~  lOOm/s  :  8.00%  Ar,  16.00%  He,  76.00%  N2 

U2  —  40  m/s  :  8.00%  Ar,  16.00%  He,  76.00%  N2  . 


Table  3.2:  Freestream  speeds  and  compositions. 
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The  local  Reynolds  numbers  for  the  shear  layers  investigated,  at  the  x  = 
0.365  m  measuring  station  (Eq.  3.6),  are  higher  than  in  the  previous  shear-layer, 
molecular-mixing  work  by  Mungal  et  al.  (1985),  which  included  an  exploratory  in¬ 
vestigation  of  inflow  effects.  In  the  facility  used  for  the  present  experiments,  it  was 
possible  to  increase  flow  speeds  (and  measuring-station  Reynolds  numbers),  such 
that  the  natural  (unperturbed)  flow  boundary  layers  could  be  expected  to  remain 
laminar,  but  at  high-enough  Reynolds  numbers  to  permit  a  transition  to  turbu¬ 
lent  boundary-layers,  if  tripped.  For  higher  Reynolds  numbers  yet,  the  naturally- 
occurring  splitter-plate  boundary  layers  would  be  transitional,  or  naturally  turbu¬ 
lent,  and  purposely  altering  inflow  conditions  by  tripping  would  not  be  expected  to 
have  as  large  an  effect. 


3.3  Results  and  discussion 

Normalized  temperature-rise  profile  measurements  for  the  natural  (untripped) 
flow  conditions  (Cases  3.1-3)  are  plotted  in  Fig.  3.1.  The  corresponding  velocity- 
profile  data  are  plotted  in  Fig.  3.2. 

To  ensure  that  the  flows  were  in  the  kinetically-fast,  high- Da  limit  (c/.  Eq.  2.5 
and  associated  discussion),  reactant  concentrations  for  the  kinetically-slower  of  the 
two  cases  (H2-rich,  <f>  =  1/8  case)  were  reduced  to  those  listed  under  Case  3.3 
(Table  3.2).  Da’s  associated  with  the  F2-rich  (<f>  =  8)  case  of  the  flip  experiment 
are  higher  than  those  of  the  H2-rich  (</  =  1/8)  case,  with  higher  concentrations 
of  the  rate-limiting  (free)  fluorine  atoms  in  the  former  case.  The  Da  ratio  for  the 
two  <j>  —  1/8  cases  is  computed  to  be  larger  than  3,  if  a  homogeneous-composition, 
constant-pressure,  “balloon”  reactor  (Dimotakis  &  Hall  1987)  is  used  to  estimate  the 
reaction-completion  times:  Tchem  =  0.45  ms  and  1.50  ms,  respectively  (c/.  Eq.  2.5). 
As  can  be  seen  in  Fig.  3.1,  the  two  <j)  =  1/8  temperature-rise  profiles,  corresponding 
to  the  already  kinetically-slower  of  the  two  ^’s,  are  within  typical  experimental 
reproducibility,  demonstrating  that  the  flows  are  in  the  high-Da  regime. 

The  values  of  the  adiabatic-flame  temperature  rise,  ATf,  for  these  runs  are  all 
in  a  regime  where  dilatation  and  other  heat-release  effects  have  been  documented  to 
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FIG.  3.1  Normalized  temperature-rise  data  with  natural  (untripped)  boundary  lay¬ 
ers.  Diamonds:  Case  3.1,  <t>  -  8.  Triangles:  Case  3.2,  <f>  =  1/8.  Asterisks: 
Case  3.3,  <f>  —  1/8,  reduced  chemical-kinetic  rate. 

be  small  (Hermanson  &  Dimotakis  1989),  once  displacement  effects  from  dilatation 
are  removed,  as  was  done  here  (see  discussion  in  Ch.  2,  above).  This  can  be  directly 
confirmed  by  comparing  the  temperature-rise  profiles,  and  1%  thicknesses,  <5t,  of 
Case  3.1  runs  with  those  of  Case  3.2  runs,  for  which  A  Tf/T0  is  substantially  different, 
by  comparing  Case  3.2  with  Case  3.3  runs,  since  both  are  in  the  fast-kinetic  regime. 
This  collapse  demonstrates  that  heat-release  effects  are  weak;  the  adiabatic  flame 
temperature  rise,  ATf,  was  lower  by  over  a  factor  of  1.5  for  the  kinetically-slower 
case  (c/.  Table  3.2).  No  significant  heat-release  effects  are  manifest  in  the  velocity 
profiles  either,  with  differences  typical  of  the  small  run-to-run  variations. 


Acoustic  forcing  is  known  to  exert  an  important  influence  on  both  large-scale 
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Fig.  3.2  Velocity  profiles  for  natural  (untripped)  flows.  Legend  as  in  Fig.  3.1. 

entrainment  ( e.g .,  Koochesfahani  &  MacKinnon  1991),  as  well  as  on  turbulent- 
fluctuation  profiles  (e.g.,  Wygnanski  et  al.  1979).  Acoustic  forcing  resonant  with 
some  local,  large-scale-passage  frequency  at  an  intermediate  downstream  location 
would  result  in  non-linear  shear-layer  growth  rates  (e.g.,  Wygnanski  et  al.  1979, 
Oster  &  Wygnanski  1982,  Ho  &  Huang  1982,  Ho  &  Huerre  1984,  Dziomba  &  Fiedler 
1985,  Roberts  1985,  and  Roberts  &  Roshko  1985). 

All  shear  layers  investigated  in  the  present  experiments  exhibit  a  very-nearly- 
linear  growth  rate,  i.e.,  6(x)  oc  x,  as  seen  in  the  flow-visualization  data  reproduced 
in  Fig.  3.3.  Such  data  represent  instantaneous  flow  visualizations,  i.e.,  are  not  time 
averaged.  They  are  subject,  however,  to  an  evolving  aspect  ratio,  b/6(x ),  where 
b  =  15.3  cm  is  the  fixed  span,  while  S(x)  oc  x,  as  the  layer  grows.  Any  three- 
dimensional  structure  of  the  flow  will  be  scaled  by  the  linearly- increasing  local  width, 
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Fig.  3.3  Composite  color-schlieren  images  of  untripped  (natural),  Case  3.2  flows. 

Note  linear  growth  and  similar  large-scale  structure,  indicating  self-similar 
flow. 


8(x),  at  its  largest  scales.  As  a  consequence,  an  apparent  evolution  in  “coherence”, 
i.e.,  an  artifact  of  the  visualization,  can  be  expected  to  make  the  flow  look  as  if  it 
is  becoming  increasingly  two-dimensional  with  increasing  downstream  distance. 


The  temperature-rise  data  in  Fig.  3.1  were  used  to  compute  the  various  quan¬ 
tities  described  above.  In  particular,  for  this  natural  (untripped)  flow: 

Case  3.1:  0  =  8,  8T/x  =  0.178,  SJ&t  =  0.277  ; 

(3.7a) 

Case  3.2:  0  =  1/8,  8T/x  =  0.177  ,  8P/8T  =  0.184  , 

where  8?  is  the  temperature-rise  profile  1%  thickness  and  <5P  is  the  product  thickness 
(Eq.  2.6).  These  numbers  are  determined  with  a  precision  for  each  flow  of  <  1%, 
with  a  run-to-run  reproducibility,  at  these  flow  conditions,  in  the  range  of  1%  —  3%. 
The  reproducibility  of  the  local  shear  layer  width,  8^/x,  between  different  runs  is 
indicative  of  this  behavior.  These  can  be  used  to  compute  the  mixed-fluid  fraction 
(Eq.  2.7), 


8m/&i  =  0.410  . 


(3.7b) 
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As  can  be  seen  in  Fig.  3.1,  the  locations  of  the  peak  temperature  rise  are  nearly- 
aligned,  despite  the  extreme  shift  in  stoichiometric-mixture  ratio  (from  <j>  =  1/8 
to  <j)  =  8).  Such  behavior  is  indicative  of  a  “non-marching”  pdf,  i.e.,  a  pdf  that 
is  nearly  independent  of  y,  corresponding  to  local  compositions  that  do  not  favor 
fluid  entrained  from  the  nearest  freestream.  These  profiles  were  recorded  in  a  high 
local  Reynolds  number  (Res  ~  2  x  105)  shear-layer  flow,  far  downstream  (x/9\  ~ 
3.3  x  103),  and  corresponding  to  a  high  value  of  the  pairing  parameter  ( P  ~  47). 
As  discussed  in  the  Introduction,  by  any  of  the  measures  employed  to  characterize 
such  flows  (cf.  Eqs.  3.1,  3.2,  and  3.3  and  3.4a-c),  these  data  were  recorded  in  an 
environment  that  must  be  regarded  as  representing  fully-developed  turbulent  shear 
flow. 


The  implied  mixed-fluid  composition  at  these  conditions  is  at  variance  with 
the  report  by  Karasso  &  Mungal  (1996).  As  discussed  by  these  authors,  one  might 
expect  some  spatial  “flip”  shift  in  the  two  temperature-rise  peaks,  i.e.,  that  the 
F2-rich  (<f>  =  8)  case  would  have  a  peak  temperature  shifted  towards  the  lean- 
reactant,  H2-bearing  freestream,  and  vice-versa  for  the  H2-rich  (<j>  =  1/8)  case,  at 
least  in  high  Reynolds  number  shear-layer  flows.  Examples  of  such  behavior  are 
found  in  the  references  cited  above  (cf.  data  and  discussion  in  Dimotakis  1991a), 
and  are  indicative  of  a  pdf  that  is  “marching”,  as  it  has  become  known,  in  which 
mixture  composition  favors  the  fluid  of  the  closer  of  the  two  freestreams  and  is 
a  weak  function  of  the  transverse  coordinate,  y.  An  absence  of  flip  shift,  similar 
to  that  seen  here,  was  observed  in  the  lower- Re,  liquid-phase,  chemically-reacting 
shear-layer  flows  of  Koochesfahani  &  Dimotakis  (1986).  Such  behavior  has  been 
attributed  to  Schmidt  number  effects,  Lagrangian- time-frame  delays  in  mixing,  or 
finite  chemical-kinetic  rate  effects  in  discussions  of  the  BBM  model  (cf.  Broadwell 
&  Mungal  1991).  Alternatively,  Rogers  &  Moser  (1992,  1994)  have  attributed  this 
to  a  higher  degree  of  two-dimensional  organization  in  the  flow.  In  their  simulations 
of  temporally-developing  shear  layers,  these  authors  found  a  correlation  between 
the  degree  of  two- dimensionality  of  conditions  used  to  initialize  the  simulation  and 
the  characteristics  of  the  pdf  that  persists  for  the  duration  of  the  simulation  time. 
In  particular,  for  initial  conditions  dominated  by  two-dimensional  disturbances, 
the  flow  structure  largely  maintained  the  two-dimensionality,  yielding  a  scalar  field 
described  by  a  non-marching  pdf.  When  initialized  by  a  more  three-dimensional 


-31- 


disturbance  field,  however,  the  flow  evolved  differently,  yielding  a  marching  pdf. 
These  temporal-flow  simulations  were  conducted  at  Reg  <  2  x  104,  i.e.,  just  past 
the  mixing  transition  (Eq.  3.1).  As  noted  by  the  authors,  the  flow  did  not  appear 
to  have  reached  a  self-similar  state  during  the  simulation  time,  at  least  for  the 
strongly-two-dimensional  case. 


y/ mm 

FIG.  3.4  Normalized  temperature-rise  data  with  tripped  high-speed  boundary  lay¬ 
ers.  Diamonds:  Case  3.1  (< j>  =  8).  Triangles:  Case  3.2  (<j>  =  1/8). 

To  explore  the  sensitivity  to  inflow  conditions,  a  0.8mm-diameter  trip  wire  was 
placed  on  the  upper  (high-speed)  face  of  the  splitter  plate,  50  mm  upstream  of  the 
splitter-plate  trailing-edge.  The  naturally-occurring  boundary-layer  is  estimated  to 
be  laminar,  with  a  momentum  thickness  of  0.095  mm  at  the  trip  location  (trailing- 
edge  values  in  Table  3.1),  estimated  using  the  method  of  Thwaites.  The  trip¬ 
wire  size  and  position  matched  those  employed  by  Mungal  et  al.  (1985),  in  the 
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FlG.  3.5  Normalized  temperature-rise  data  for  tripped  and  untripped  boundary  lay¬ 
ers  at  <f>  =  8.  Both  experiments  with  Case  3.1  reactant  concentrations. 
Diamonds:  untripped  high-speed  boundary  layer,  triangles:  tripped  high¬ 
speed  boundary  layer. 


previous  (subsonic)  version  of  this  facility.  White  (1974)  recommends  a  Reynolds 
number  based  on  wire  diameter  in  excess  of  0.85  x  103  to  ensure  transition  (in  a 
receptive  boundary  layer).  The  trip  wire  Reynolds  number  in  this  flow  is  estimated 
to  be  4.6  x  103.  The  resulting  temperature-rise  and  streamwise- velocity  profile 
data,  corresponding  to  those  plotted  in  Figs.  3.1  and  3.2,  are  plotted  in  Figs.  3.4 
and  3.7,  respectively,  with  schlieren  data  reproduced  in  Fig.  3.8.  In  addition,  a 
direct  comparison  of  the  normalized  temperature-rise  profiles  in  the  untripped  and 
tripped  flows  is  found  in  Figs.  3.5  and  3.6,  for  the  F2-rich  (<f>  =  8)  and  H2-rich 
((f>  =  1/8)  stoichiometries,  respectively. 
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FIG.  3.6  Normalized  temperature-rise  data  for  tripped  and  untripped  splitter-plate 
boundary  layers  at  <j>  =  1/8.  Both  experiments  with  Case  3.2  reactant 
concentrations.  Legend  as  in  Fig.  3.5. 


The  temperature-rise  profiles  in  Fig.  3.4  are  measured  0.415  m,  or,  approxi¬ 
mately,  500  trip-wire  diameters,  downstream  of  the  trip  location;  in  the  “far-field”  of 
this  perturbation.  The  untripped  flow,  at  the  same  location,  satisfies  all  previously- 
proposed  criteria  for  fully- developed  flow.  In  the  case  of  the  tripped  boundary-layer 
flow,  even  if  the  boundary  layer  (momentum  thickness)  was  assumed  to  increase  to 
1/3  the  trip-wire  diameter,  the  same  criteria  would  still  be  satisfied. 


The  shear-layer  growth  and  mixed-fluid  quantities,  for  the  flows  with  the  high¬ 
speed  boundary-layer  tripped,  were  also  derived  from  these  temperature-rise  profile 
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Fig.  3.7  Velocity  profiles  for  the  flows  with  tripped  high-speed  side  boundary  layers. 
Legend  as  Fig.  3.4. 

data, 

Case  3.1:  <£  =  8  ,  =  0.140  ,  ^p/^T  =  0.293  ; 

(3.8a) 

Case  3.2:  <£  =  1/8,  d^/x  =  0.138  ,  <5P/<5T  =  0.211  , 

yielding  a  mixed-fluid  fraction  of, 

6m/6T  =  0.448  .  (3.8b) 

Comparing  with  the  untripped  boundary-layer  data  (Eq.  3.7),  tripping  the  high¬ 
speed  boundary  layer  can  be  seen  to  have  resulted  in  a  growth  rate,  S(x)  j x.  that  is 
lower  by  21%  ±  1.5%.  This  must  be  regarded  as  a  substantial  reduction  considering 
that  this  flow  should  have  attained  quasi-asymptotic  behavior,  i.e.,  one  in  which 


FIG.  3.8  Composite  color-schlieren  images  of  tripped  (0.8mm  wire  on  high-speed 
side),  Case  3.2  flows.  Note  linear  growth  and  geometrically-similar  large- 
scale  structures. 

inflow/initial  conditions  are  expected  to  exert  a  minimal  influence,  if  any.  This  re¬ 
duction,  however,  is  in  qualitative  accord  with  previous  investigations,  at  lower  Reg, 
of  two-stream  (as  opposed  to  half-jet)  shear  layers,  as  discussed  in  the  Introduction. 

Significantly,  the  temperature-rise  data  for  the  tripped  high-speed  boundary- 
layer  exhibit  a  “flip”  shift,  indicating  that  a  transition  to  a  “tilted”  molecularly- 
mixed-fluid  pdf  has  occurred,  i.e.,  the  mixed-fluid  composition  is  now  a  function 
of  the  transverse  coordinate,  y.  If  this  is  to  be  expected  of  turbulence  in  a  “more- 
developed”  state,  as  suggested  by  Karasso  &  Mungal  (1996),  it  cannot  be  attributed 
to  an  increase  in  Reynolds  number.  The  21%  reduction  in  local  thickness,  8(x), 
induced  by  tripping  the  high-speed  side  boundary  layer,  decreases  the  local  Reynolds 
number,  Reg,  by  the  same  amount.  The  data  are  also  at  variance  with  the  proposed 
Karasso  &  Mungal  criteria  for  pdf  evolution  and  behavior,  in  terms  of  the  change 
in  x j 6\ .  The  tripped  flow  has  a  slightly-thicker  high-speed  boundary  layer,  i.e., 
a  higher  6\,  and,  therefore,  a  lower  scaled  downstream  distance,  x/d\,  as  well  as 
pairing  parameter,  P.  Changes  in  this  direction  should  have  had  the  opposite 
effect,  i.e.,  taken  a  “non-marching”  pdf,  characteristic  of  “undeveloped”  flow,  to  a 
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Fig.  3.9  Normalized  temperature-rise  data  for  tripped  and  untripped  splitter-plate 
boundary  layers  at  <j)  —  8.  Transverse  coordinate  normalized  by  <5t-  Dia¬ 
monds:  untripped  high-speed  side,  triangles:  tripped  high-speed  side. 

“marching”  pdf,  characteristic  of  “developed”  flow,  according  to  these  authors. 

Tripping  reduced  entrainment  and  shear-layer  growth  far  downstream  of  the 
change  in  inflow  conditions.  Interestingly,  it  also  produced  more  vigorous  stirring 
and  effective  molecular  mixing.  The  latter  can  be  ascertained  by  plotting  the  data, 
scaling  the  transverse  coordinate  by  the  local  width,  i.e.,  as  in  Figs.  3.9  and 

3.10.  It  is  in  these  units  that  the  computation  of  the  mixed-fluid  fraction,  Sm/S-£,  is 
performed  (Eq.  2.7).  The  increase  in  chemical  product  and  molecular  mixing  within 
the  shear  layer,  in  each  case,  is  clear.  As  can  be  seen  by  comparing  the  product 
fractions,  6p/8 t,  for  the  two  cases  (Eqs.  3.7  vs.  3.8),  the  increase  is  approximately 
6%  in  the  <f>  =  8  flow  and  16%  in  the  <f>  =  1/8  flow. 
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Fig.  3.10  Normalized  temperature-rise  data  for  tripped  and  untripped  boundary 
layers  at  <f>  =  1/8.  Transverse  coordinate  normalized  by  <$t-  Legend  as 
Fig.  3.9. 

Comparing  the  schlieren  image  data  in  Fig.  3.3  and  Fig.  3.8  reveals  that  tripping 
the  high-speed-side  boundary  layer  has  not  only  changed  the  growth  rate,  but  also 
the  large-scale  structure,  far  downstream.  For  the  tripped  flow,  a  decrease  in  two- 
dimensional  organization  of  the  flow  is  evident,  with  spanwise  optical  integration 
obscuring  local  large-scale-structures  in  the  upstream  half  of  the  flow.  The  inferred 
“marching”  scalar  pdf  behavior  for  this  flow  is  observed  along  with  this  decrease 
in  two-dimensional  organization,  as  also  noted  by  Rogers  &  Moser  (1994)  in  their 
numerical  simulations.  In  the  context  of  large-scale-entrainment  models  for  this  flow 
(Dimotakis  1986,  1991b),  a  change  in  the  entrained-fluid  mole  ratio,  Fn,  is  expected. 
While  this  was  not  measured,  one  can  infer  the  mixed-fluid  molar  composition  ratio, 
estimated  with  Eq.  2.8.  The  two  ratios,  En  and  £n>  need  not  have  the  same  value. 
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FlG.  3.11  Velocity  profiles  for  both  untripped  and  tripped  flows.  Transverse  coor¬ 
dinate  normalized  by  1%  temperature-rise  thickness. 


The  former  measures  the  entrained-fluid  ratio,  at  the  largest  scales  of  the  flow, 
whereas  the  latter  measures  the  mixed-fluid  ratio.  Using  the  tripped  and  untripped 
shear  layer  temperature-rise  measurements,  we  find  (Eq.  2.8), 

(£n)untripped  -  L52  >  V$ ■  (Stripped  “  L39  >  (3-9) 

i.e.,  a  shift  to  a  lower  mixed-fluid  ratio  (closer  to  symmetric).  Interestingly,  the 
increase  in  the  large-scale-structure  spacing-to-position  ratio  induced  by  tripping, 
evident  in  the  schlieren  data,  suggests  that  the  change  in  En  is  towards  larger  values, 
i.e.,  in  the  opposite  direction. 

The  velocity  profiles  for  the  untripped  and  tripped  flows  examined  thus  far  axe 
plotted  in  Fig.  3.11,  vs.  the  normalized  transverse  coordinate,  y/8 t,  for  comparison. 
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Fig.  3.12  Normalized  temperature- rise  data  with  various  low-speed  boundary  layer 
trips.  All  flows  with  Case  3.2  reactant  concentrations.  Diamonds:  Un¬ 
tripped  flow.  Triangles:  0.8  mm  trip  wire.  Asterisks:  1.6  mm  trip  wire. 

This  scaling  can  be  seen  to  produce  a  substantial  collapse  of  the  mean,  streamwise 
velocity  profile  at  the  measuring  station,  despite  the  differences  between  tripped 
and  untripped  flows  discussed  above.  This  collapse  also  demonstrates  that  the 
transverse  length  scale,  as  estimated  on  the  basis  of  mixed-fluid/temperature-rise 
data,  also  scales  the  velocity  field,  i.e.,  that  <5u(x)  ~  <5T(x),  where  6u(x)  would  be 
some  appropriate  velocity-profile-based  transverse  width.  Finally,  the  collapse  in  the 
scaled  profiles  in  Fig.  3.11  indicates  that  the  mean  velocity  is  largely  insensitive  to 
the  processes  responsible  for  the  changes  in  entrainment  and  mixing.  Accordingly,  it 
would  be  surprising  if  shear-layer-entrainment  and  -mixing  models  based  on  mean- 
velocity  profiles,  e.g.,  of  eddy- viscosity-type,  could  accommodate,  let  alone  predict, 
the  observed  mixing  dependence  on  inflow  conditions. 
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FlG.  3.13  Normalized  temperature-rise  data  recorded  with  various  boundary  layer 
trips.  Both  flows  employ  Case  3.2  reactant  concentrations,  and  0.8  mm- 
diameter  trip  wires.  Diamonds:  Both  sides  tripped,  triangles:  High-speed 
side  tripped. 

Additional  experiments  were  performed  to  explore  the  effects  that  other  near¬ 
field  perturbations  might  have,  at  these  nominal  conditions.  Guided  by  the  results 
of  Lang  (1985),  in  which  a  large  low-speed-side  momentum  defect  (wake)  was  found 
to  significantly  increase  the  shear-layer  growth  rate  at  a  freestream  speed  ratio  of 
r  ~  0.5,  two  different  trip  wires  were  employed  on  the  low-speed  side  of  the  splitter 
plate,  with  diameters  of  0.8  mm  and  1.6  mm,  also  at  x  =  —50  mm.  The  natural 
low-speed  side  boundary  layer  is  estimated  to  be  laminar  and  to  have  a  momentum 
thickness  of  0.15  mm  at  this  trip  location.  As  seen  in  Fig.  3.12,  these  low-speed-side 
perturbations  have  little,  if  any,  influence  on  the  far-field  growth  of  the  layer,  in 
contrast  to  the  significant  changes  induced  by  tripping  the  high-speed  splitter  plate 
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boundary  layer,  and  only  a  weak  effect  on  mixing. 

A  possible  reason  for  the  lack  of  low-speed-side  boundary-layer  trip  influence 
is  that  the  naturally-occurring  (untripped)  boundary  layer  is  already  turbulent. 
The  low-speed  side  splitter  plate  has  concave  curvature  over  a  portion  of  its  length 
(Fig.  2.3),  and  may  be  susceptible  to  Taylor-Gortler  vortex  formation.  At  the  con¬ 
ditions  of  these  experiments,  however,  the  Taylor  number  is  estimated  to  be,  ap¬ 
proximately,  4,  significantly  less  than  the  stability  limit  of  7  proposed  by  Liepmann 
(1945;  cf.  also  discussion  in  Schlichting  1968).  This  indicates  that  the  untripped 
low-speed  boundary  layer  did  not  support  Taylor-Gortler  vortices  and  was  likely 
laminar. 

As  a  further  exploration,  both  splitter-plate  boundary-layers  were  tripped  to  see 
if  additional  effects  could  be  detected  (cf.  Fig.  3.13).  In  this  case,  the  temperature- 
rise  profile  is  nearly  indistinguishable  from  that  of  the  flow  in  which  only  the  high¬ 
speed  side  layer  was  tripped. 

The  lack  of  shear-layer  sensitivity  to  low-speed  side  perturbations,  relative  to 
that  seen  for  high-speed  perturbations,  is  likely  due  to  the  stability  characteristics 
of  the  near-field  flow.  In  particular,  linear-stability  analyses  of  wake-like  profiles 
exhibit  a  “resonance”  phenomenon,  similar  to  that  found  in  quantum  mechanics 
(Koochesfahani  &  Freiler  1989).  In  their  study  of  density-ratio  effects  on  wake¬ 
like  flows,  these  authors  found  substantial  increases  in  amplification  of  the  “wake” 
mode,  relative  to  the  usually-amplified  “shear-layer”  mode,  for  particular  choices 
of  the  mean  profile  corresponding  to  such  a  resonance  condition.  A  similar  effect 
may  be  present  in  this  s  ~  1  flow,  where  the  relative  amplification  of  various  modes 
might  be  substantially  altered  by  a  change  in  the  initial  velocity  profile  (tripping). 
In  particular,  the  profile  for  the  high-speed-tripped  flow  may  be  closer,  in  a  state- 
space  sense,  to  such  a  resonance  than  the  low-speed-tripped  profile.  If  such  a  mode 
were  amplified  to  the  non-linear  regime,  the  flow  might  find  itself  in  a  different 
“basin  of  attraction”,  a  behavior  observed  in  low- dimensional  dynamical  systems, 
e.g .,  Wiggins  (1990). 

Table  3.3  summarizes  the  various  computed  quantities  for  each  of  the  natural, 
or  tripped,  inflow  boundary  layer  configurations. 
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Mixture,  perturbation 

6? /mm 

6p/mm 

<W^T 

Res/ 10 

Case  3.1,  untripped: 

65.0 

18.0 

0.277 

2.3 

Case  3.1,  0.8  mm  HSS  trip: 

51.1 

15.0 

0.293 

1.8 

Case  3.2,  untripped: 

64.6 

11.9 

0.184 

2.3 

Case  3.2,  0.8  mm  HSS  trip: 

50.5 

10.7 

0.211 

1.8 

Case  3.2,  0.8  mm  LSS  trip: 

63.7 

12.0 

0.189 

2.2 

Case  3.2,  1.6  mm  LSS  trip: 

63.7 

12.4 

0.195 

2.2 

Case  3.2,  0.8  mm  trips  (both): 

51.3 

10.6 

0.206 

1.8 

Case  3.3,  untripped: 

67.4 

12.3 

0.181 

2.4 

(  /^T  )untripped 

~  0.410 

(^m/^T  )tripped 

~  0.448 

(^n)untripped  —  1.52 

(^n)  tripped  —  1.39 

Table  3.3:  Chemically-reacting  flow  data  for  various  inflow  conditions. 

3.4  Conclusions 

In  the  shear-layer  flows  investigated,  the  measuring  station  is  well  downstream 
of  the  minimum  recommended  distance,  as  scaled  by  the  initial  boundary-layer (s) 
momentum  thickness,  the  values  of  the  “pairing  parameter”  (Eq.  3.4)  are  well  in 
excess  of  previously-proposed  minimum  values,  and  the  local  Reynolds  numbers  are 
well  beyond  the  “mixing  transition”,  i.e.,  Res  »  104,  as  indicated  for  asymptotic, 
self-similar  behavior,  mixing,  and  chemical-product  formation.  In  this  flow  environ¬ 
ment,  the  shear-layer  growth  rate,  6t(x)/x,  the  mixed-fluid  mole  fraction,  6m/6T, 
as  well  as  the  mixed-fluid  composition  ratio,  £n,  were  all  found  to  depend  on  inflow 
conditions.  The  present  results  corroborate  previously-documented  findings  regard¬ 
ing  inflow  (boundary-layer)  effects  on  shear-layer  growth.  They  also  demonstrate, 
however,  that  inflow  conditions  influence  small-scale  (molecular)  mixing  in  shear 
layers,  far  downstream  and  at  high  values  of  the  local  Reynolds  number.  Interest¬ 
ingly,  scaling  the  velocity  profiles  downstream  by  shear-layer  thickness  leaves  no 
other  effect  of  inflow  conditions  on  these  profiles,  concealing  the  conspicuous  effects 
on  (molecular)  mixing. 


The  data  presented  here,  as  well  as  in  previous  investigations,  indicate  that 
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while  a  quasi-self-similar  state  can  be  achieved,  at  high  Reynolds  numbers,  in  a 
relatively  short  fraction  of  the  downstream  distance,  this  state  is  a  function  of 
the  inflow  conditions.  The  flow  at  high  Reynolds  numbers  and  away  from  inflow 
boundaries  is  not  found  to  be  a  unique  function  of  local-flow  parameters,  such  as 
freestream- velocity  and  -density  ratios,  local  Reynolds  number,  etc.  No  support 
for  the  existence  of  a  universal,  local-parameter-defined  flow  could  be  gleaned  from 
these  investigations.  Indeed,  the  apparent  discrepancies  between  the  new  results 
and  previous  reports,  as  well  as  inconsistencies  between  numerous  previous  reports, 
are  probably  attributable  to  interpretations  that  have  implicitly  assumed  that  an 
inflow-independent,  quasi-asymptotic  state  exists. 
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CHAPTER  4 


Compressibility  and  dynamics 


This  chapter  details  results  from  a  set  of  non-reacting  experiments  in  the 
transonic-  and  supersonic-flow  regime.  Interactions  between  pressure  waves,  both 
stationary  and  traveling  in  the  laboratory  frame,  and  turbulence  are  present  in 
all  flows  investigated.  Freestream  Mach  numbers  were  varied  to  produce  both 
bi-supersonic  and  bi-subsonic  flows,  to  explore  what  effect,  if  any,  the  presence 
of  (lab-frame)  subsonic  regions  might  have  on  the  flow.  Among  other  findings, 
turbulent-structure  convection  speeds  were  observed  to  violate  an  empirical  mode- 
selection  rule  that  was  based  on  the  presence  of  a  subsonic  freestream.  In  addition, 
a  sensitivity  of  the  global  flowfield  to  various  boundary  conditions  was  observed, 
due  primarily  to  transonic-flow  sensitivity  to  imposed  area  changes. 
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4.1  Introduction 

Compressibility  is  an  important  parameter  of  turbulent  shear-layer  flows.  As 
described  in  Ch.  1,  experimental  evidence  demonstrates  that  an  increase  in  velocity 
difference  between  the  two  freestreams,  as  scaled  by  a  representative  speed  of  sound, 
is  correlated  with  a  decrease  in  shear-layer  growth  rate,  when  normalized  by  the 
corresponding  incompressible  value  for  the  same  freestream-velocity  and  -density 
ratios  ( e.g .,  Bogdanoff  1983,  Chinzei  et  al.  1986,  Papamoschou  &  Roshko  1988, 
Samimy  &  Elliott  1990,  Goebel  &  Dutton  1991,  Clemens  &  Mungal  1992a,  and 
Hall  et  al.  1993). 

Compressible-flow  data  also  indicate  that  the  two-dimensional,  coherent  struc¬ 
tures  observed  in  subsonic  shear  layers  do  not  persist  as  compressibility  is  in¬ 
creased.  In  particular,  Clemens  &  Mungal  (1992a,  1995)  noted  that  their  images, 
obtained  using  Mie  scattering  from  an  alcohol  condensate  as  well  as  nitric-oxide 
planar  induced  fluorescence,  indicate  strong  (large-scale)  three-dimensionality,  for 
Mc  >  0.62.  These  experiments,  however,  covered  a  streamwise  extent  in  the 

early  portion  of  the  shear-layer  flow,  with  guidewalls  separated  by  a  vertical  dis¬ 
tance,  h/8(x )  »  1,  i.e.,  large  when  scaled  by  the  local  shear-layer  width.  In  this 
regime,  the  shear  layer  can  be  expected  to  approach  unbounded-flow  behavior.  Data 
recorded  at  Mc  ~  0.5  (Hall  et  al.  1993,  Figs.  lOa-c)  indicate  a  smooth  transition 
from  a  flow  close  to  the  splitter  plate,  where  h/8(x)  >  1,  in  which  spanwise  (two- 
dimensional)  organization  is  absent,  to  one  further  downstream,  where  h/8(x)  ~  1, 
where  a  two-dimensional,  large-scale  organization  emerges.  Similar  behavior  is  ex¬ 
hibited  in  the  downstream  part  of  the  higher-Mc  spanwise-integrated  image  data  in 
Clemens  &  Mungal  (1995,  Fig.  31).  The  degree  of  spanwise  coherence  can  be  further 
enhanced  downstream  via  the  interaction  with  a  strong,  2-D  wave  system,  especially 
when  a  resonant-wave  condition  is  encountered  between  large-scale- structure  spac- 
ings  and  streamwise  round-trip  distances  of  wall-reflected  waves  (Tam  &  Hu  1989; 
and  Hall  et  al.  1993,  Figs.  11).  The  experimental  data  are,  however,  subject  to  an 
evolving  aspect  ratio,  b/8(x),  as  the  shear-layer  grows  (see  discussion  in  Ch.3),  an 
effect  that  would  favor  visualization  of  an  apparent  enhanced  coherence  at  large  x, 
i.e.,  small  b/8(x). 

Linear-stability  analyses  of  unbounded  flows  show  suppression  of  the  growth 
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o£  two-dimensional,  spanwise-coherent  modes  and  a  relative  increase  in  the  growth 
of  three-dimensional,  oblique  waves  for  Mc  >  0.6.  In  the  presence  of  (reflecting) 
walls,  however,  where  most  of  these  flows  have  necessarily  been  experimentally 
investigated,  the  amplification  of  two-dimensional  disturbances  can  be  dominant 
( e.g .,  Tam  &:  Hu  1989,  and  Zhuang  et  al.  1990). 

As  described  in  Ch.  1,  one  of  the  characteristics  of  moderate-to-high  compress¬ 
ibility  shear-layer  flows  are  supersonically-convecting  features  that  produce  Mach- 
wave  radiation  in  the  freest  reams.  One  such  convected-wave  system  was  captured 
by  Hall  et  al.  (1993,  see  Fig.  5  and  related  discussion),  in  the  subsonic  freestream  of 
a  high-compressibility  shear  layer.  They  estimated  (by  measuring  the  wave-system 
Mach  angles)  a  convective  Mach  number  of  the  supersonic  disturbance  with  respect 
to  the  low-speed  freestream,  in  the  range, 

2-1  <  Mc2  <  2.6  ,  at  Mc  ~  0.96 ,  (4.1) 

indicating  the  presence  of  a  fast-mode  disturbance  in  the  flow.  As  noted  by  those 
authors,  the  waves  are  essentially  acoustic  in  nature  (near-specular  reflection  from 
the  lower  guidewall)  and  do  not  appear  to  be  associated  with  a  two-dimensional, 
spanwise-coherent,  large-  ( 6 -)  scale  structure.  Despite  the  relative  lack  of  flow  struc¬ 
ture  observed  in  direct- visualization  experiments  at  similar  conditions,  these  waves 
are  evidently  produced  by  a  large-scale,  long-lived  (temporally-coherent)  distur¬ 
bance,  as  they  are  seen  to  emanate  with  a  characteristic  spacing  scaled  by  the  local 
outer  scale,  6(x ),  and  propagate  for  long  distances  in  a  near-linear  manner. 

It  has  been  recently  proposed  that  these  waves  are  a  result  of  bounding-wall 
influences  (e.g.,  Umemura  et  al.  1996  and  Smits  &  Dussuage  1996).  To  explore  this 
possibility,  we  refer  to  the  original  visualizations  of  Hall  et  al.  (1993),  where  the 
waves  are  seen  to  be  generated  immediately  following  the  splitter-plate  tip,  where 
the  flow  is  essentially  unbounded,  i.e.,  h/6(x)  >  1.  The  radiation  is  continually 
emitted  by  the  layer  to  the  test-section  exit,  where  the  flow  is  likely  to  be  partially 
influenced  by  the  walls,  since  h/6(x )  ~  1  there.  The  characteristic  Mach  angle 
of  the  emitted  waves  remains  approximately  constant  over  this  large  range  of  h/6, 
indicating  that  the  convection  speed  of  the  wave-generating  turbulent  features  is 
independent  of  the  bounding-wall  geometry,  at  least  in  that  flow. 
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We  also  note  that  such  a  wave  system  does  not  preclude  the  mode  coexistence 
observed  by  Oertel  (1979),  since  the  fast  mode  would  likely  be  the  only  source  of 
turbulence-generated  acoustic  radiation  in  this  flow  (Dimotakis  1991b).  Preliminary 
data  reported  in  Dimotakis  &  Leonard  (1991),  recorded  at  high-compressibility 
conditions  (Mc  c*  1.5),  appear  to  indicate  the  presence  of  two  distinct  modes,  each 
radiating  into  a  separate  freestream. 

As  indicated  by  Eq.4.1,  observations  of  turbulent-structure  convection  speeds, 
or  equivalently,  convective  Mach  numbers,  are  at  variance  with  (central-mode) 
isentropic-recovery  estimates,  t.g.,  Eq.  1.5.  Instead,  features  are  found  to  con- 
vect  with  a  speed  close  to  one  of  the  freestream  speeds.  Proposals  by  both  Pa- 
pamoschou  (1991)  and  Dimotakis  (1991a, b)  indicated  that  the  ability  of  shear-layer 
structures  to  support  convected  shock  waves  could  be  responsible  for  this  asymme¬ 
try.  The  proposal  of  Dimotakis  included  a  model  calculation  of  convective  Mach 
numbers  that  took  total-pressure  (shock)  losses  into  account,  and  was  found  to  be 
in  quantitative  accord  with  existing  data. 

Papamoschou  (1989)  proposed  that  this  convection-speed  bifurcation  is  gov¬ 
erned  by  the  presence  (or  absence)  of  a  laboratory-frame  subsonic  freestream,  i.e., 

Supersonic/subsonic  flow  <=>  fast-mode  convection  speed,  i.e.,  Mc\  <C  Mc 2, 


and 


Bi-supersonic  flow  <=>  slow-mode  convection  speed,  i.e.,  Mc  1  »  Mc 2. 

Subsequent  supporting  data,  e.g.,  Hall  (1991),  led  to  the  coining  of  this  proposal 
as  a  “stream-selection  rule”  by  Dimotakis  (1991b),  who,  however,  noted  its  diffi¬ 
culty  with  respect  to  Galilean  invariance.  Obviously,  a  (purely)  convective-frame 
description  of  the  flow  cannot  capture  a  criterion  of  this  nature. 

We  also  note,  with  regard  to  downstream-to-upstream  communication,  e.g., 
from  the  exit  region  of  the  test-section  to  the  splitter-plate  tip,  that  an  elliptic 
(subsonic)  freestream  provides  an  acoustic-feedback  path,  possibly  coupling  the 
near-field  of  shear-layer  development  to  disturbances  generated  far  downstream. 
The  precise  effect  of  such  feedback  is  not  clear.  However,  as  discussed  by  Dimotakis 
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&  Brown  (1976),  it  potentially  exercises  a  considerable  influence  in  incompressible 
flows.  This  communication  path  is  absent  in  a  bi-supersonic  configuration,  effec¬ 
tively  decoupling  upstream  (inflow)  conditions  from  the  facility  outflow  conditions. 


4.2  Experiments 

Experimental  flow  conditions  for  this  part  of  the  study  are  listed  in  Table 
4.1.  Evident  in  many  of  the  bi-supersonic  flow  visualizations,  e.g.,  Fig.  4.1,  are 
spatially-varying  freestream  conditions  close  to  the  test-section  inlet,  at  small  x. 
These  variations  are  due  to  differences  between  freestream  static  pressures  and  flow 
directions  at  the  splitter- plate  tip,  in  addition  to  (discontinuous)  changes  in  flow  di¬ 
rection  induced  by  shear-/mixing-layer  displacement  thicknesses.  In  a  bi-supersonic 
configuration,  the  flow  can  only  avail  itself  of  these  differences  through  a  system  of 
shock  and  expansion  waves.  Measurements  presented  in  Table  4.2,  e.g.,  Sv-is/x,  are 
recorded  in  the  far  field,  once  the  two  freestream  flows  have  “equilibrated”,  however, 
the  mechanisms  and  dynamics  of  this  near-field  equilibration  are  often  noteworthy. 
As  before,  guidewalls  were  adjusted  to  provide  non-accelerating  freestream  flows  in 
the  downstream  measurement  region. 


Gas  1 

Gas  2 

7i 

72 

Mle 

M2e 

Case  4.1 

n2 

n2 

1.40 

1.40 

1.14 

0.65 

Case  4.2 

n2 

n2 

1.40 

1.40 

1.42 

1.13 

Case  4.3 

Ar 

n2 

1.67 

1.40 

1.47 

1.22 

Case  4.4 

He 

n2 

1.67 

1.40 

1.39 

1.12 

Case  4.5 

He 

(Na/C^)1 

1.67 

1.30 

0.9 

0.5 

1In  equal  (molar)  proportions. 

Table  4.1:  Freestream  gas  properties  and  nominal  flow  conditions. 


The  term  “mixing  layer”  is  used  to  describe  some  flows,  since  low-shear  (r  ~ 
1)  layers  are  likely  influenced,  at  least  in  early  regions  of  development,  by  the 
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splitter-plate  wake  (Zhuang  &:  Dimotakis  1995).  Additionally,  when  used  to  describe 
pressure  waves,  the  terms  “stationary”  and  “traveling”  refer  to  a  laboratory-fixed 
reference  frame. 

As  indicated  in  Eq.  1.6,  presentation  of  compressibility  effects  on  growth  rate  re¬ 
quires  normalization  by  an  estimate  of  the  corresponding  (same  r,  s )  incompressible- 
flow  growth  rate.  This  normalization  is  provided  by  Eq.  1.2,  appropriate  for  the 
spatially-developing  layers  investigated  here. 


4.3  Results  and  discussion 


4.3.1  Low-compressibility  flows 

To  explore  bi-supersonic  flow,  i.e.,  flow  with  two  laboratory-frame  supersonic 
freestreams,  a  set  of  rims  was  initially  performed  at  low-compressibility,  i.e.,  Mc  ~ 
0.1.  Significant  dynamical  differences  from  bi-subsonic  flow  are  inherent  to  this 
regime,  namely  the  absence  of  downstream-to-upstream  feedback  and  the  presence 
of  refractive  interactions  between  stationary  pressure  waves  and  turbulence.  Such 
interactions  are  disallowed,  i.e.,  are  purely  reflective,  in  supersonic/subsonic  flow. 

To  investigate  Case  4.1,  the  facility  was  configured  to  produce  a  subsonic  top 
freestream  and  a  supersonic  bottom  freestream  at  the  test-section  inlet  plane.  Both 
gases  were  pure  N2.  As  seen  in  the  schlieren  image  in  Fig.  4.1,  this  configuration 
results  in  a  flow  different  from  all  others  presented  here;  the  lower  freestream  is  the 
high-speed  stream.  Evident  in  this  visualization  is  the  interaction  of  the  mixing 
layer  with  a  fan  from  the  5°  expansion  corner  at  the  (upstream)  test-section  inlet. 
The  subsequent  turning  of  the  instantaneous  dividing  streamline  initially  causes  a 
reduction  in  the  cross-sectional  area  of  the  top,  subsonic  freestream  and  a  corre¬ 
sponding  increase  in  the  cross-sectional  area  of  the  bottom,  supersonic  freestream, 
effectively  producing  an  aerodynamic  Laval  nozzle. 
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Fig.  4.1  Schlieren  (vertical  knife-edge)  visualization  of  Case  4.1  mixing-layer  flow. 

As  predicted  by  quasi-one-dimensional  gasdynamics,  the  Mach  numbers  and, 
consequently,  the  velocities  of  both  freestreams  are  increased  in  this  situation.  There 
is,  however,  no  gain  in  total  momentum  flux,  as  the  increase  in  dynamic  pressure 
is  offset  by  a  corresponding  decrease  in  static  pressure,  as  required  by  momentum 
considerations.  Although  the  top  freestream  is  subsonic  at  the  inlet  plane,  it  accel¬ 
erates  to  a  supersonic  condition,  as  evidenced  by  the  near-vertical  wave  transmitted 
into  the  top  freestream.  This  top-freestream  acceleration  continues  for  several  wave 
reflections,  until  the  freestream  static  pressures  and  flow  directions  equilibrate,  to 
form  a  bi-supersonic  flow.#  The  turbulence  does  not  appear  to  be  influenced  by 
the  transition  from  a  subsonic  to  supersonic  freestream  region,  other  than  (minor) 
growth-rate  variations  attributable  to  local  freestream-velocity  variations. 

An  interesting  shockwave-interaction  phenomenon  is  seen  in  the  top  freestream, 
downstream  from  the  transition  to  supersonic  flow.  A  reflection  of  a  compression 
wave  from  the  top  guidewall  appears  to  be  a  Mach  reflection.  However,  as  discussed 
by  various  authors,  e.g .,  Hornung  (1986),  such  a  configuration  is  disallowed  by  the 


This  hybrid  freestream  condition  was  originally  dubbed  a  “bi-sonic”  flow,  until  Dr.  Julian 
Tishkoff  pointed  out  that  the  experiment  had  nothing  to  do  with  the  large  animals  that  once 
roamed  the  plains. 
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classic  (3-shock)  von  Neumann  theory,  at  this  (local)  freestream  Mach  number  of 
1.13.  The  resolution  of  this  “von  Neumann  paradox”  is  a  matter  of  current  study 
and  debate.  Proposed  explanations  concentrate  on  the  triple-point  region,  with 
modifications  to  the  three-shock  theory  including  slipstream  displacement  effects 
(Kobayashi  et  al.  1995),  finite-curvature  effects  (Colella  &  Henderson  1990,  and 
Brio  &  Hunter  1992),  and,  perhaps  the  most  likely  explanation,  the  addition  of 
a  fourth,  centered-expansion  wave  (Guderley  1947,  see  discussion  in  Sandemann 
1997).  The  wave  reflection  observed  in  Fig.  4.1  is  an  example  of  such  an  interaction 
in  steady  flow,  in  contrast  to  pseudo-steady  (self-similar  in  mixed  spatio-temporal 
variables)  flows  generated  in  shock-tube  experiments. 

A  second  set  of  low-compressibility  mixing-layer  flows  are  visualized  in  Figs.  4.2. 
These  flows  are  distinguished  from  Case  4.1  above,  in  that  both  freestreams  are 
supersonic  at  the  nozzle-exit  plane.  Evident  in  the  visualization  of  Case  4.2a 
(Fig.  4.2a),  however,  is  a  near-normal,  stationary,  compression  shock  in  the  bot¬ 
tom  freestream,  located  slightly  downstream  of  the  nozzle-exit  plane.  An  oblique 
Mach  wave  incident  on  the  normal-wave  front  is  seen  to  terminate  there,  as  de¬ 
manded  by  post-shock  subsonic  (wave-frame)  flow.  The  flow  is  then  accelerated 
by  the  effective  freestream-area  change  produced  by  the  combination  of  mixing- 
and  boundary-layer  displacement  thicknesses.  Notably,  the  mixing  layer  presents  a 
negative  displacement  thickness  to  the  bottom  freestream  flow,  while  the  boundary 
layer  presents  a  positive  displacement.  The  net  effect  is  a  transition  back  to  su¬ 
personic  flow  a  short  distance  downstream  of  the  normal  wave,  as  is  evident  from 
the  refractive  interaction  of  a  (stationary)  Mach  wave  with  the  mixing  layer,  at 
that  downstream  location.  This  post-shock  acceleration  is  a  product  of  the  same 
mechanism,  i.  e.,  an  aerodynamic  Laval  nozzle,  by  which  the  top  freestream  of  Case 
4.1  transitioned  to  a  bi-supersonic  configuration. 

This  response  to  local  displacement  thicknesses  is  a  consequence  of  the  large  ef¬ 
fect  that  small  cross-sectional-area  changes  can  have  in  transonic  flow,  and  is  largely 
due  to  the  sensitivity  of  the  quasi-one-dimensional  Mach  number-area  relationship 
in  this  regime  (e.g.,  Figure  5.2  of  Liepmann  &  Roshko  1957).  An  excellent  example 
of  this  sensitivity  is  contained  in  the  transonic-airfoil  experiments  of  Liepmann  (see 
sequence  245-248  in  Van  Dyke  1982),  where  small  differences  between  displacement 


FlG.  4.2a  Schlieren  (oblique  knife  edge)  visualization  of  bi-supersonic  Case  4.2a 
mixing-layer  flow. 


FlG.  4.2b  Schlieren  (oblique  knife  edge)  visualization  of  bi-supersonic  Case  4.2b 
mixing-layer  flow. 

thicknesses  of  laminar  and  turbulent  boundary  layers  produce  large  changes  in  the 
global  airfoil  flowfield. 

In  Case  4.2a,  there  is  a  transition  from  supersonic  to  subsonic,  and  eventually 
back  to  supersonic,  flow.  This  generates  a  subsonic  (lab-frame)  freestream  region 
in  a  fashion  similar  to  the  Case  4.1  flow,  but  here  the  subsonic  region  is  separated 
from  the  inflow  conditions.  Once  again  there  is  no  qualitative  change  in  turbulence 
properties  associated  with  transitions  between  supersonic  and  subsonic  lab-frame 
freestream  flow. 

A  perturbation  of  Case  4.2a  flow  conditions  results  in  Case  4.2b.  The  latter  is 
produced  with  a  slightly-lower  (5%)  high-speed  freestream  stagnation  pressure.  This 
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perturbation  of  the  inlet-plane  “boundary  condition”  results  in  a  flow  configuration 
shown  in  Fig.  4.2b,  where  the  normal  shock  observed  previously  is  now  “unfolded” 
into  a  Mach  reflection  from  the  bottom  guidewall.  Visualizations  recorded  in  a 
sequence  of  decreasing  stagnation  pressures  exhibit  a  gradual  shortening  of  the 
Mach  stem,  until  regular  reflection  is  attained.  As  illustrated  by  Case  4.2  flows, 
it  is  possible  to  manipulate  the  strength  of  these  waves  by  adjusting  upstream 
boundary  conditions,  e.g.,  stagnation  pressures.  Because  of  the  discontinuous  jump 
in  displacement  thicknesses  at  the  splitter-plate  tip,  however,  such  a  wave  system 
will  always  be  generated,  even  as  an  infinitesimally-weak  N-wave. 

These  flow  configurations,  i.e.,  two  moderately-supersonic,  acoustically-similar 
gas  streams,  provide  a  unique  diagnostic  capability.  As  noted  above,  bi-supersonic 
flow  permits  refraction  of  stationary  waves  by/through  the  mixing  layer.  In  the  flows 
discussed  above,  turbulence-generated  index-of-refraction  (density)  gradients  axe 
relatively  weak  in  comparison  to  wave-generated  gradients.  As  observed  in  Fig.  4.2a, 
for  example,  this  discrepancy  in  the  two  types  of  gradients  provides  instantaneous 
visualizations  of  wave-turbulence  interactions.  Given  that  many  of  these  waves 
axe  well-appxoximated  as  infinitesimal  (acoustic)  distuxbances,  such  data  serve  as 
passive  measurements  of  the  instantaneous  Mach  number  field.  Interestingly,  this 
field  appears  to  be  nearly  discontinuous  in  many  cases,  i.e.,  waves  axe  “kinked”. 
This  suggests  an  instantaneous  Mach  number  field  quite  different  from  the  smoothly- 
varying  mean  field,  a  behavior  also  observed  in  incompressible  flow  (Dimotakis  et 
al.  1981). 

In  the  case  where  these  interactions  are  of  finite  strength  (non- acoustic),  the 
potential  for  baroclinic-vorticity  generation  exists.  Such  vorticity  would  be  aligned 
with  the  mixing-layer  vorticity,  either  augmenting  or  diminishing  it,  depending  on 
the  density  gradient.  This  interaction  geometry  is  purely  two-dimensional,  in  con¬ 
trast  to  the  experiment  of  Clemens  &  Mungal  (1992b)  in  which  enhanced  (spanwise- 
integrated)  layer  growth  was  observed  as  a  result  of  a  three-dimensional  shock¬ 
wave/shear-layer  interaction.  In  such  interactions,  the  baroclinically-generated  vor¬ 
ticity  is  not  parallel  to  that  of  the  original  mixing  layer,  and  alters  (direction  and 
magnitude)  the  vorticity  field,  creating  a  substantially-different  flowfield. 

In  N2/N2  flows  above,  the  potential  for  vorticity  generation  was  not  large,  as 
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Fig.  4.3  Color  schlieren  photograph  of  Case  4.3  mixing-layer  flow. 


freestream-density  gradients  are  close  to  unity.  To  explore  the  effects  on  a  layer 
that  supports  a  significant  density  gradient,  a  flow  similar  to  Case  4.2  was  inves¬ 
tigated,  but  with  a  high-speed  freestream  of  Ar,  instead  of  N2.  The  low-speed 
freestream  in  Case  4.3  was  still  N2,  resulting  in  freestream  conditions  of  s  ~  0.54 
and  r  ~  0.98  (at  the  shock-impingement  location),  conditions  that  could  result  in 
significant  (fractional)  augmentation  of  mixing-layer  vorticity.  However,  as  seen 
in  the  color  schlieren  photograph  in  Fig.  4.3,  there  is  no  discernible  shock-induced 
change  in  the  mixing-layer  growth  rate.  It  should  be  noted  that  the  unperturbed 
vorticity  field  in  this  flow  is  probably  wake-like  in  nature,  as  the  relative  shear  is 
low.  This  finding  is  in  accord  with  (two-dimensional)  shock- wave/mixing-layer  in¬ 
teraction  experiments  of  other  investigators,  e.g.,  Shau  et  al.  (1993),  who  concluded 
that  a  similar  interaction  had  little  net  effect  on  large  scales  in  the  flow. 


4.3.2  Moderate-compressibility  flows 

In  contrast  to  the  low-compressibility  flows  discussed  above,  Case  4.4  and  4.5 
flows  were  produced  with  helium  high-speed  freestreams.  In  Case  4.4,  for  example, 
the  relatively-high  speed  of  sound  of  helium  yields  a  higher  compressibility  than 
that  of  Case  4.2,  at  the  same  nominal  freestream  Mach  numbers  (Table  4.1).  This 
feature  was  originally  exploited,  in  this  facility,  by  Hall  (1991)  and  Hall  et  al.  (1993). 
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Two  distinct  pressure-wave  systems  are  evident  in  Case  4.4  (Figs.  4.4ab)  flow 
visualizations,  a  stationary-wave  system  similar  to  that  found  in  low-compressibility 
flows,  and  a  traveling-wave  system  in  the  low-speed  freestream,  similar  to  that  ob¬ 
served  by  Hall  et  al.  (1993).  As  in  Case  4. 1-4.3  flows,  the  stationary-wave/shear- 
layer  interaction  does  not  significantly  modify  the  growth  rate,  although  this  in¬ 
teraction  is  estimated  to  produce  vorticity  of  relatively  small  magnitude  and  of 
opposite  sign  to  the  unperturbed  shear-layer  vorticity.  In  addition,  this  wave  sys¬ 
tem  does  not  appear  to  be  resonant  with  any  two-dimensional  structure  in  the  flow; 
the  enhanced  two-dimensional  structure  observed  in  the  downstream  region  by  Hall 
et  al.  (1993)  is  not  observed  here  (Fig.  4.4b).  We  note,  however,  that  the  values  of 
h/8  in  this  flow  are  less  than  those  for  enhanced-structure  flow,  and  Case  4.4  may 
approximate  unbounded-flow  conditions. 


FlG.  4.4a  Schlieren  (oblique  knife  edge)  photograph  of  the  upstream  region  of  bi- 
supersonic  Case  4.4  shear-layer  flow. 

As  in  the  Hall  et  al.  (1993)  analysis  (Eq.  4.1),  a  convective  Mach  number  with 
respect  to  the  low-speed  freestream  is  inferred  from  the  image  data  as, 

Case  4.4  :  MC2  —  1.70  ±  0.10  =£•  Mci  ~  0.29  ±  0.017.  (4-2) 

This  is  in  general  accord  with  the  previously-observed  asymmetry,  i.e.,  Mc i  ^  Mc 2, 
as  shown  in  Fig.  4.5.  Only  measurements  recorded  in  two-dimensional  shear-layer 
flows,  he.,  no  data  from  annular-shear-layer/turbulent-jet  flows,  are  presented  here. 
Recent  evidence  indicates  that  geometrical  differences  are  likely  to  be  important 
(Freund  et  al.  1997). 


FlG. 4.4b  Schlieren  (oblique  knife  edge)  photograph  of  downstream  region  of  bi- 
supersonic  Case  4.4  flow.  Same  nozzle-exit  conditions  as  those  recorded 
in  Fig.  4.4a. 


FlG.  4.4c  Schlieren  (oblique  knife  edge)  photograph  of  off-design  (now  bi-subsonic) 
Case  4.4  flow. 

Although  in  general  agreement  with  the  asymmetry,  Case  4.4  exhibits  a  fast¬ 
mode  convection  speed  in  bi-supersonic  flow;  a  result  at  variance  with  the  stream- 
selection  rule  originally  proposed  by  Papamoschou  (1989).  An  exception  to  this 
empirical  rule  is  also  in  the  original  Papamoschou  data  (1989,  Case  N31N17),  how¬ 
ever,  at  an  ambiguously-low  level  of  compressibility. 

This  measurement  provides  direct  evidence,  at  a  moderate-to-high  level  of  com¬ 
pressibility  (Mc  ~  0.7),  that  the  (empirically-based)  “stream-selection  rule”,  a  cri¬ 
terion  based  solely  on  laboratory-frame  conditions,  does  not  provide  an  accurate 
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FIG. 4.5  Observed  turbulent-feature  convection  speeds.  Isentropic-recovery  model 
predicts  a  line  with  near-unity  slope,  i.e.,  Mc\  ~  Mc i- 

description  of  the  convection-speed  behavior.  This  is  manifest  in  Fig.  4.6,  where  the 
data  have  been  split  into  groups  according  to  whether  the  freestream  flows  are  super¬ 
sonic,  or  not.  As  with  the  low-compressibility  results  discussed  earlier,  the  presence 
of  (here,  the  lack  of)  a  subsonic  laboratory-frame  freestream  region  does  not  govern 
large-scale  dynamics  in  turbulent  shear/mixing  layers.  Several  possible  alternative 
criteria,  e.g.,  an  analog  of  the  regular-/Mach-reflection  “information  condition”  of 
Homung  et  al  (1979),  were  also  tested  against  the  data  with  no  success. 

The  sensitivity  to  imposed  freestream  boundary  conditions  exhibited  by  the 
low-compressibility  flows  was  also  observed  here.  In  particular,  the  flow  conditions 
shown  in  Fig.  4.4c  were  recorded  during  the  multiple-run  iterative  procedure  re¬ 
quired  to  produce  (on-design)  bi-supersonic  Case  4.4  flow  conditions.  This  flow  was 
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Fig.  4.6  Stream-selection  data.  Case  4.4:  bi-supersonic  flow  on  vertical  (fast-mode) 
branch. 

produced  with  the  same  (nominal)  nozzle-/contraction-plenum  pressures,  but  at  a 
non-diverging  guidewall  angle,  i.e.,  0°  as  opposed  to  1.2°  divergence  for  on-design 
Case  4.4,  as  is  evident  from  the  comparison  of  Fig.  4.4c  and  Figs.  4.4ab,  respectively. 
The  flow  is  almost  completely  subsonic  in  this  off-design  configuration,  having  been 
shocked  by  a  three- wave  structure  near  the  test-section  inlet.  Interestingly,  a  trav¬ 
eling  wave  system  is  again  present  in  the  low-speed  freestream,  similar  to  that 
observed  in  bi-supersonic  Case  4.4  flow. 

The  possibility  for  supersonically-convecting  turbulent  structures  in  bi-subsonic 
flow  is  further  addressed  by  Case  4.5  conditions.  As  seen  in  Fig.  4.7,  there  is  indeed 
a  traveling  (oblique)  pressure-wave  system  in  the  bottom  (low-speed)  freestream, 


FlG.  4.7  Schlieren  (vertical  knife  edge)  photograph  of  Case  4.5  bi-subsonic  shear- 
layer  flow. 

indicative  of  a  disturbance  with  convective  Mach  numbers, 

Case  4.5  :  Mc2  c*  1.40  ±  0.25  =►  Mcl  ~  0.20  ±  0.036.  (4.3) 

There  is  also  an  upstream-traveling,  i.e.,  right-to-left,  set  of  near-normal  pres¬ 
sure  waves,  similar  to  those  observed  by  Hall  (1991)  in  the  same  facility.  These 
disturbances  are  the  signature  of  a  downstream-to-upstream  feedback  mechanism 
which  exploits  subsonic  flow.  The  probable  origin  of  this  acoustic  feedback  is  the 
large-scale-structure  communication  described  by  Dimotakis  &  Brown  (1976).  As 
noted  by  these  authors,  this  feedback  is  present  in  any  subsonic-flow  facility,  cou¬ 
pling  outflow  (downstream)  conditions  to  those  in  the  inflow  region. 

Normalized  growth  rates  for  these  two  flows,  as  well  as  those  of  two  low- 
compressibility  flows  (Cases  4.1  and  4.2),  along  with  previous  data,  are  presented 
in  Table  4.2  and  Fig.  4.8.  Only  data  recorded  in  two-dimensional  shear-layer  exper¬ 
iments  are  plotted,  for  the  same  geometrical  reasons  detailed  for  convection-speed 
data.  Absolute  growth  rates  of  Case  4.1  and  4.2  flows  are  extremely  small,  and  as 
a  consequence,  have  significant  fractional  uncertainties  associated  with  them,  be¬ 
cause  of  near-unity  values  of  the  velocity  ratio,  r.  Case  3.3  is  excluded  from  this 
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Fig.  4.8  Normalized  compressible  shear- layer  growth  rates  as  parameterized  by  the 
total  convective  Mach  number,  Mc  (Eq.  1.5).  Normalization  provided  by 
Eq.  1.2,  with  Cg's  as  suggested  (when  documented)  by  the  various  investi¬ 
gators. 

presentation,  as  it  is  more  accurately  classified  as  a  wake  than  a  shear-layer.  The 
data  are  seen  to  be  in  good  overall  agreement  with  those  of  previous  investigations, 
suggesting  that  the  presence  of  subsonic/supersonic  freestream-fiow  regions  does 
not  have  a  significant  effect  on  growth  rate. 
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M1 

m2 

r 

s 

Mc 

60/x 

*^vis/ 

vis /^( 

Case  4.1: 

1.31 

1.06 

0.85 

0.91 

0.097 

0.0275 

0.028 

1.02 

Case  4.2: 

1.41 

1.09 

0.82 

0.94 

0.15 

0.034 

0.035 

1.04 

Case  4.3:1 

1.44 

1.26 

0.98 

0.54 

- 

- 

- 

Case  4.4: 

1.49 

1.27 

0.34 

5.30 

0.70 

0.072 

0.22 

0.33 

Case  4.5: 

0.9 

0.5 

0.21 

5.80 

0.52 

0.150 

0.33 

0.46 

1  Conditions  at  shock-impingement  location. 

Table  4.2:  Far-held  (equilibrated)  flow  properties. 

4.4  Conclusions 

Non-reacting  shear- /mixing-layer  flows,  in  the  transonic  regime,  have  been  in¬ 
vestigated.  Interactions  between  pressure  waves  and  turbulence  are  manifest  in 
schlieren  visualizations.  Three  flows,  with  transonic  laboratory-frame  freestreams, 
but  at  low  levels  of  mixing-layer  compressibility  (Mc  ~  0.1),  were  found  to  be  sen¬ 
sitive  to  imposed  boundary  conditions.  Both  bi-supersonic  and  bi-subsonic  flows, 
at  Mc  ~  0.7  and  Mc  ~  0.5,  respectively,  are  found  to  exhibit  turbulent-structure 
convection  speeds  that  violate  a  previously-proposed  rule  for  mode  selection.  The 
moderate-compressibility,  bi-supersonic  flowfield  was  also  observed  to  be  sensitive 
to  imposed  boundary  conditions,  further  demonstrating  the  potential  for  global- 
flowfield  manipulation,  with  significant  control  authority,  in  this  transonic-flow 
regime. 

Normalized  growth  rates  measured  in  these  flows  are  in  qualitative  accord  with 
existing  data;  regardless  of  whether  the  flow  is  bi-supersonic,  bi-subsonic,  or  a 
hybrid-flow  condition.  Along  with  the  violation  of  the  stream-selection  rule,  a  crite¬ 
rion  based  on  lab-frame  Mach  numbers,  this  suggests  that  the  presence,  or  absence, 
of  laboratory-frame  subsonic  freestream  regions  does  not  significantly  affect  free- 
shear  turbulence. 
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CHAPTER  5 


Compressibility  and  mixing 


This  chapter  describes  experimental  results  regarding  compressibility  effects 
on  turbulent  shear-layer  mixing.  Two  chemically-reacting  “flip”  experiments  were 
performed,  while  holding  freestream-velocity  and  -density  ratios,  Reynolds  num¬ 
ber,  and  reactant  concentrations  nominally  fixed.  An  increase  in  compressibility, 
while  holding  these  other  parameters  fixed,  is  associated  with  a  decrease  in  growth 
rate,  in  accord  with  previous  investigations.  The  increase  in  compressibility  is  also 
associated  with  an  increase  in  the  peak  (normalized)  temperature  rise,  and  to  a 
lesser  degree,  the  mixed-fluid  fraction.  This  change  in  small-scale  mixing  is  solely 
attributable  to  compressibility. 
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5.1  Introduction 

It  is  often  stated  that  the  effect  of  compressibility  on  turbulent  shear-layer  flows 
is  to  suppress  mixing ,  e.g.,  Martens  et  al.  (1994).  Such  an  inference  is  likely  made 
on  the  basis  of  data  concerning  the  suppression  of  the  (large-scale)  growth  rate  by 
compressibility.  However,  as  discussed  by  various  authors,  e.g.,  Dimotakis  (1986), 
high-Reynolds-number  mixing  processes  can  be  conceptually  broken  into  three  sep¬ 
arate  stages:  large-scale  entrainment,  intermediate-scale  stirring,  and  small-scale 
diffusion/mixing,  with  each  stage  important  to  the  overall  molecular  mixing  of  the 
two,  initially-separated,  freestream  constituents. 

This  sequence  can  be  expressed  in  terms  of  the  quantities  that  are  accessible 
to  the  measurement  techniques  employed  here,  i.e., 


The  first  term  represents  the  first  stage  of  mixing,  while  the  second  represents  the 
process  at  all  of  the  remaining  spatial  scales,  down  to  the  diffusion,  or  Batchelor, 
scale.  The  behavior  of  6/x  as  a  function  of  compressibility  is  relatively  well  known; 
the  behavior  of  Sm/S  (here  <5m/<5x)  is  a  matter  of  some  controversy. 

The  partition  in  Eq.  5.1  indicates  that  we  should  expect  a  net  decrease  in  the 
amount  of  molecularly-mixed  fluid,  at  least  at  moderate-to-high  compressibility 
conditions.  This  expectation  stems  from  two  experimental  facts  (Dimotakis  1991a): 
the  large-scale  growth  rate  is  approximately  1/2  of  the  incompressible  value  by  Mc  ~ 
0.5  ( e.g.,  Ch.  4),  and  the  mixed-fluid  fraction  is  approximately  1/2  at  incompressible- 
flow  conditions  {e.g.,  Ch.  3). 

There  have  been  several  recent  reports  on  molecular  mixing  in  compressible 
shear-layer  flows,  i.e.,  Hall  et  al.  (1991),  Clemens  Sz  Paul  (1995),  Island  (1997),  and 
Freund  et  al.  (1997),  with  an  almost  equal  number  of  differing  conclusions.  The 
techniques  employed  in  these  studies  are  threefold:  chemically-reacting  flip  experi¬ 
ments  (this  work,  and  Hall  et  al.  1991),  cold-chemistry  (fluorescence-quenching)  flip 
experiments  (Clemens  &  Paul  1995,  Island  1997),  and  direct  numerical  simulation 
(Freund  et  al.  1997).  This  variety  of  techniques,  in  addition  to  differing  measures 
of  6m/6  employed  in  each  case,  renders  direct  comparison  difficult. 
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Hall  et  al.  (1991)  concluded  that  the  amount  of  molecularly-mixed  fluid  de¬ 
creased  with  increasing  Mc,  Clemens  &  Paul  (1995)  concluded  that  it  increased, 
Island  (1997)  concluded  no  change  or  a  slight  increase,  and  Freund  et  al.  (1997) 
observed  an  increase.  The  three  sets  of  experimental  data  are,  however,  subject 
to  several  caveats.  In  particular,  changes  in  Mc  are  accompanied  by  associated 
changes  in  velocity  ratio,  density  ratio,  and/or  Reynolds  number.  As  a  result,  in¬ 
terpretations  of  the  data  are  often  ambiguous,  with  (strong)  assumptions  required. 
A  typical  example  of  such  an  ambiguity  is  found  in  the  Island  (1997)  interpreta¬ 
tion,  where  it  is  unclear  whether  observed  increases  in  8m/6  are  attributable  to 
compressibility  effects,  or  Reynolds-number  effects,  or  even  possible  large-scale  ef¬ 
fects,  e.g.,  streamwise  pressure  gradient,  and/or  velocity /density  ratio.  The  Hall  et 
al.  (1991)  and  Clemens  &  Paul  (1995)  data  are  also  subject  to  similar  ambiguities. 
The  numerical  simulations  of  Freund  et  al.  (1997)  are  free  of  such  difficulties,  as  all 
parameters  except  Mach  number  were  held  nominally  fixed.  The  Reynolds  number 
of  the  simulated  flow  is  however,  very  low  because  of  computational  constraints. 
In  particular,  the  simulated-flow  Reynolds  numbers  of  Re  ~  2.5  x  103  (based  on 
jet  radius)  make  direct  comparison  with  experimentally-generated  flows  difficult,  as 
the  experiments  are  conducted  at  Reynolds  numbers  well  in  excess  of  the  mixing 
transition. 

Collectively,  this  evidence  suggests  that  if  compressibility  effects,  as  measured 
by  the  convective  Mach  number,  Mc,  are  to  be  assessed  with  any  confidence,  we 
must  hold  as  many  of  these  parameters  as  constant  as  possible.  Such  a  set  of 
shear-layer  flip  experiments  is  described  below. 


5.2  Experiments 

These  experiments  were  performed  at  two  nominal  sets  of  freestream  condi¬ 
tions,  a  low-compressibility  flow  with  Mc  —  0.25  (Cases  5.Xa)  and  a  moderate- 
compressibility  flow  with  Mc  cz  0.47  (Cases  5.Xb). 

Freestream-velocity  and  -density  ratios  were  chosen  to  facilitate  a  comparison 
with  the  existing  body  of  chemically-reacting  data  (see  discussion  in  Ch.  3,  and 
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5.1a  -  Fluorine-rich  (<j>  =  8,  ATf  ~  171 K): 

lh  ~  280  m/s  :  8.00%  Ar,  15.07%  He,  0.93%  H2,  0.14%  NO,  75.86%  N2 

U2  ^  110  m/s  :  8.00%  Ar,  8.00%  He,  8.00%  F2,  76.00%  N2  . 

5.2a  -  Hydrogen-rich  {<j>  =  1/8,  ATf  ~  311 K): 

Ux  ~  280  m/s  :  8.00%  Ar,  14.00%  H2,  0.88%  NO,  77.12%  N2 

U2  —  110 m/s  :  14.20%  Ar,  8.00%  He,  1.80%  F2,  76.00%  N2  . 

5.3a  -  Hydrogen-rich,  increased  concentration  (4>  =  1/8,  ATf  ~  355  K): 

TJ\  ~  280  m/s  :  8.00%  Ar,  16.00%  H2,  1.00%  NO,  75.00%  N2 

U2  ~  llOm/s  :  13.94%  Ar,  8.00%  He,  2.06%  F2,  76.00%  N2  . 

5.4a  -  Non-reacting: 

Ux  ~  280  m/s  :  8.00%  Ar,  16.00%  He,  76.00%  N2 
U2  ~  110  m/s  :  16.00%  Ar,  8.00%  He,  76.00%  N2  . 

Table  5.1a:  Low-compressibility  flow:  freestream  speeds  and  compositions, 
references  therein),  i.e., 

r  =  ^  ~  0.4 ,  and  s  =  —  —  1.0 .  (5-2) 

Ux  Pi 

Low-compressibility  data  were  recorded  in  flows  with  freestream  speeds, 

Ux  —  280 m/s  and  U2  ~  110 m/s  =>  AU  —  170 m/s,  (5.3a) 

and  moderate-compressibility  data  in  flows  with, 

Ux  ^  525 m/s  and  U2  ~  210 m/s  AU  ~  315 m/s.  (5.3b) 

In  both  flows,  freestream  compositions  documented  in  Tables  5.1  were  tailored  to 
offset  the  various  density  changes  associated  with  compressible-flow  freestreams, 
i.e.,  a  freestream  molecular-mass  ratio  of  unity  does  not  imply  s  =  1.  In  addition, 
the  flip-experiment  mixture  ratio,  <f>o,  was  held  constant  at  <f> o  =  8. 
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5.1b  -  Fluorine-rich  (<f>  =  8,  ATf  ~  179  K): 

U\  ~  525  m/s  :  41.07%  He,  0.93%  H2,  0.14%  NO,  57.86%  N2 

U2  ~  210 m/s  :  20.00%  Ar,  14.00%  He,  8.00%  F2,  58.00%  N2  . 

5.2b  -  Hydrogen-rich  ((f)  =  1/8,  ATf  ~  329  K): 

Ux  ~  525  m/s  :  28.00%  He,  14.00%  H2 ,  0.88%  NO,  57.12%  N2 

U2  m  210  m/s  :  26.20%  Ar,  14.00%  He,  1.80%  F2,  58.00%  N2  . 


5.3b  -  Hydrogen-rich,  increased  concentration  (<f>  =  1/8,  ATf  ~  375  K): 

Ux  ~  525  m/s  :  26.00%  He,  16.00%  H2,  1.00%  NO,  57.00%  N2 

U2  ~  210 m/s  :  25.94%  Ar,  14.00%  He,  2.06%  F2,  58.00%  N2  . 


5.4b  -  Non-reacting: 

Ux  ^  525  m/s  :  42.00%  He,  58.00%  N2 

U2  ~  210  m/s  :  28.00%  Ar,  14.00%  He,  58.00%  N2  . 


Table  5.1b:  Moderate-compressibility  flow:  freestream  speeds  and  compositions. 


With  reference  to  Eq.  5.3,  it  might  appear  that  the  Reynolds  number  of  the 
moderate-compressibility  flow  would  be  significantly  higher  than  that  of  the  low- 
compressibility  flow.  As  discussed  above,  such  an  increase  has  been  a  shortcoming 
of  previous  investigations.  The  present  experiment  was  designed  to  exploit  the 
compressibility-induced  reduction  in  growth  rate  (for  example,  in  Fig.  4.8),  to  pro¬ 
duce  minimal  changes  in  Reynolds  number  as  compressibility  was  changed.  There 
are  also  relatively  small  changes  in  mixture  density  and  viscosity  associated  with 
lower  (static)  temperatures  in  the  moderate-compressibility  flow. 


Profiles  of  the  temperature  rise  from  chemical  reaction  were  recorded  at  a 
downstream  station  of  x  =  0.365  m,  as  before.  Splitter-plate  boundary  layers  in 
these  flows  are  thinner  (higher  Reg)  than  for  flows  described  in  Ch.  3,  which  were 
themselves  seen  to  comfortably  satisfy  all  indicated  criteria  for  fully-developed,  self¬ 
similar  flow.  The  resulting  shear  layers  have  a  local  Reynolds  number  range  at  this 
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measuring  station  of, 


7.6  x  105  <  Res  <  7.9  x  105  , 


(5.4a) 


for  Case  5.Xa  flows  and, 


9.8  x  105  <  Re6  <  10  x  105  ,  (5.4b) 

for  Case  5.Xb  flows.  The  slight  (1/10-decade)  change  in  Res  is  due  to  imposed 
requirements  discussed  above,  i.e .,  constant  r,  s,  and  <j> o. 

The  increase  in  velocity  does,  however,  result  in  a  substantial  increase  of  the 
total  convective  Mach  number,  from, 

Mc  ~  0.25  ,  (5.5a) 


for  Case  5.Xa  flows,  to, 


Mc  ~  0.47  , 


(5.5b) 


for  Case  5.Xb  flows.  As  a  consequence  of  this  restriction  of  the  parameter  space, 
any  changes  in  mixed-fluid  properties  can  be  attributed  solely  to  the  change  in 
compressibility. 


5.3  Results  and  discussion 

Normalized  temperature-rise  profile  measurements  for  low-compressibility  flow 
conditions  (Cases  5.1a-5.3a)  axe  plotted  in  Fig.  5.1. 

A  procedure  similar  to  that  detailed  in  Ch.  3  was  employed  to  ensure  that  the 
flows  were  in  the  kinetically-fast  regime.  In  this  case,  however,  reactant  concen¬ 
trations  for  the  kinetically- slower  H2-rich  flows  were  increased,  in  contrast  to  the 
reduction  employed  in  Ch.  3  flows.  The  Da  ratios  for  the  <f>  =  1/8  cases,  in  both 
low-  and  mo  derate- compressibility  flows,  are  computed  to  be  larger  than  1.5.  Once 
again,  Damkohler  numbers  associated  with  the  F2-rich  (<f>  =  8)  part  of  the  flip 
experiment  are  considerably  higher  than  those  of  the  H2-rich  (<j>  =  1/8)  part. 


y/mm 


FlG.  5.1  Normalized  temperature-rise  data  for  the  low-compressibility  flow  (Mc  ~ 
0.25).  Diamonds:  Case  5.1a,  (f>  =  8.  Triangles:  Case  5.2a,  (j)  =  1/8. 
Asterisks:  Case  5.3a,  <f>  =  1/8,  increased  chemical-kinetic  rate. 


Fig.  5.2  Composite  color-schlieren  image  of  Case  5.2a  flow.  Note  linear  growth  and 
similar  large-scale  structure. 


The  lower-concentration  compositions  were  used  to  estimate  the  mixed-fluid 
quantities,  to  alleviate  concerns  regarding  the  effects  of  heat  release.  In  particular, 
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at  these  high-speed  conditions,  mixing  times,  rmjx,  are  substantially  reduced  from 
those  in  Ch.  3  flows,  requiring  a  commensurate  reduction  in  chemical  times,  rchenu 
to  remain  in  the  fast-kinetic  regime.  This  is  accomplished  by  increasing  reactant 
concentrations  in  the  H2-rich  (<f)  =  1/8)  case.  There  are,  of  course,  associated 
increases  in  adiabatic  flame  temperatures  and  consequent  levels  of  heat  release.  As 
seen  in  Fig.  5.1,  there  may  be  slight  heat-release  effects  at  work;  the  temperature 
profile  recorded  for  (higher-concentration)  Case  5.3a  concentrations  is  reduced  in 
both  width  and  height  relative  to  that  recorded  for  Case  5.2a,  indicative  of  slight 
reductions  in  both  shear-layer  growth  and  product  fraction,  respectively.  These 
reductions  are  in  accord  with  those  observed  by  Hermanson  &  Dimotakis  (1989), 
at  similar  heat-release  levels. 


Fig.  5.4  Composite  color-schlieren  image  of  Case  5.2b  flow. 


The  temperature-rise  data  in  Fig.  5.1  were  used  to  compute  the  various  quan¬ 
tities  described  above.  In  particular,  for  this  low- compressibility  flow, 

Case  5.1a:  <j)  =  8  ,  6^/x=  0.174  ,  <5p/<$t  =  0.287  ; 

(5.6a) 

Case  5.2a:  (f>  =  1/8  ,  St/x  =  0.164  ,  £p/#t  =  0.229  , 

which  in  turn  yield  a  mixed-fluid  (mole)  fraction, 

6m/ST  =  0.459  ,  (5.6b) 


and  a  mixed-fluid  composition  ratio, 


Sn  = 


1.25  . 


(5.6b) 
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Fig.  5.3  Normalized  temperature-rise  data  for  moderate-compressibility  flow  (Mc  ~ 
0.47).  Same  scale  as  in  Fig.  5.1,  above,  to  allow  direct  comparison  of  the 
profiles.  Diamonds:  Case  5.1b,  <f>  =  8.  Triangles:  Case  5.2b,  <f>  =  1/8. 
Asterisks:  Case  5.3b,  <f>  =  1/8,  increased  chemical-kinetic  rate. 


Normalized  temperature-rise  measurements  for  moderate-compressibility  flow 
conditions  (Cases  5. lb-5. 3b)  are  plotted  in  Fig.  5.3.  Several  differences  relative  to 
low-compressibility  flow  (Fig.  5.1)  are  noted.  In  particular,  the  shear-layer  is  seen 
to  be  thinner,  i.e.,  6/x  has  been  reduced,  at  this  higher  level  of  compressibility, 
in  quantitative  agreement  with  previous  investigations  (Fig.  4.8).  For  the  H2-rich 
(<f>  =  1/8)  parts  of  the  flip  experiment,  the  maximum  normalized  temperature  rises 
are  comparable  in  both  low-  and  moderate-compressibility  flows.  However,  the 
maximum  attained  in  the  F2-rich  (<f>  =  8)  part  is  significantly  (25%)  higher  than  in 
low-compressibility  flow,  reaching  nearly  70%  of  the  adiabatic  flame  temperature. 
This  behavior  was  also  observed  in  moderately-compressible  flows  with  freestream- 
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FlG. 5.5  Normalized  temperature-rise  data  for  F2-rich  (<f>  =  8)  flows.  Diamonds: 

low-compressibility  (Mc  ~  0.25)  flow,  triangles:  moderate-compressibility 
flow  (Mc  ~  0.47). 

diluent  mixtures  of  helium  and  argon  (Bond  1996,  private  communication). 

As  in  Eqs.  5.6  above, 


Case  5.1b: 

<f>  =  8  ,  8t/x  =  0.103  , 

Sp/6T  =  0.303  ; 

(5.7a) 

Case  5.2b: 

(f>  =  1/8  ,  S'l/x  =  0.107  , 

8v/8r  =  0.221  , 

which  in  turn  yield 

a  mixed-fluid  fraction, 

8m/8r  =  0.466  , 

(5.7b) 

and  a  mixed-fluid  composition  ratio, 


(5.7b) 


=  1.37  . 
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These  data,  despite  the  significant  increase  in  normalized  temperature-rise  at 
4>  =  8,  show  a  only  a  modest  increase  in  chemical  product  fraction,  and  conse¬ 
quently,  mixed-fluid  fraction,  associated  with  the  increase  in  Mc.  This  finding 
reflects  a  change  in  the  shape  of  the  temperature-rise  profile,  and  consequently, 
the  mixed-fluid  structure.  Figure  5.5  provides  a  direct  comparison  of  the  F2-rich 
temperature-rise  profiles,  with  the  transverse  coordinate  normalized  by  Sy  (see  dis¬ 
cussion  in  Ch.  3.  Changes  are  also  observed  in  the  instantaneous  flow  visualizations 
in  Figs.  5.2  and  5.4,  where  qualitative  differences  in  visualized  structure  axe  evident 
A  modified  flow  structure  is  also  indicated  provided  by  the  change  in  the  mixed- 
fluid  composition  ratio,  which  is  seen  to  increase  for  moderate-compressibility  flow, 
favoring  the  high-speed  fluid  to  a  greater  extent  than  in  low-compressibility  flow. 


5.4  Conclusions 

Effects  of  compressibility  on  molecular  mixing  were  investigated  with  the  use 
of  chemically-reacting  flip  experiments  at  two  different  freestream-flow  conditions. 
These  conditions  are  characterized  by  nominally-fixed  local-flow  Reynolds  numbers, 
freestream-velocity  and  -density  ratios,  and  reactant  compositions,  but,  a  change  in 
convective  Mach  number  from  Mc  ~  0.25  to  Mc  ~  0.47. 

This  (isolated)  increase  in  flow  compressibility  allows  the  observed  changes  in 
shear-layer  growth  and  mixing  to  be  attributed  solely  to  the  increase  in  Mc.  A  25% 
increase  in  the  peak  normalized  temperature  rise  was  recorded  in  the  F2-rich  part 
of  the  flip  experiment.  However,  compressibility  has  also  altered  the  molecularly- 
mixed  fluid  structure,  as  indicated  by  changes  in  temperature-rise  profile  shapes,  as 
well  as  in  the  mixed-fluid  composition  ratios.  As  a  result,  the  associated  increase 
in  mixed-fluid  fraction  is  small. 
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CHAPTER  6 


An  alternative  compressibility  parameter 


A  new  shear-layer  flow  compressibility  parameter  is  proposed  as  an  alterna¬ 
tive  to  the  (total)  convective  Mach  number,  Mc.  This  parameter  is  derived  from 
consideration  of  compressibility  as  an  energy  transformation  (kinetic  to  thermal, 
and  vice  versa )  mechanism,  and  differs  from  Mc  at  non-unity  freestream-density 
(speed-of-sound)  ratios.  Experimentally-observed  growth  rates  are  well  represented 
by  this  new  parameter. 
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6.1  Introduction 

Compressibility  effects  on  shear-layer  growth  are  often  presented  in  terms  of 
the  total  convective  Mach  number,  in  the  form  indicated  by  Eq.  1.6.  This  scaling 
captures  the  overall  reduction  in  growth  rate  with  increasing  Mc,  as  discussed  ear¬ 
lier.  Figure  4.8  is  reproduced  below  as  Fig.  6.1a,  for  reference.  A  fair  degree  of 
“scatter”  is  evident  in  the  data.  However,  as  noted  by  Clemens  &  Mungal  (1995), 
among  others,  the  scaling  exhibited  by  many  individual  sets  of  experimental  data, 
e.g .,  Chinzei  et  al.  (1986),  is  reasonably  consistent.  This  behavior  is  reminiscent  of 
that  documented  in  Ch.  3,  i.e.,  an  influence  of  inflow  (or  other  facility-dependent) 
conditions  that  result  in  the  range  of  values  exhibited  by  Eq.  1.3.  Such  a  variance  is 
not  typically  addressed  in  reports  of  compressible  shear-layer  growth  rates,  where 
a  value  for  C$  may  (often)  be  chosen  to  match  experimental  data,  e.g.,  of  Brown 
&  Roshko  (1974).  An  additional  complicating  difficulty  is  presented  by  the  use  of 
different  measures  of  the  local  thickness,  8,  i.e.,  pitot  thickness,  visual  thickness, 
etc.,  that  comprise  the  data  shown  in  Fig.  6.1a,  and  the  ill-defined  relationships 
between  various  measures  (Papamoschou  1986). 

In  an  attempt  to  reconcile  experiments  performed  in  various  facilities,  the  con¬ 
stant  C(  is  re-evaluated  for  each  individual  experimental  data  set.  This  procedure, 
and  the  resulting  values  for  Cs,  are  documented  in  Appendix  B.  As  seen  in  Fig.  6.1b, 
this  renormalization  provides  a  substantial  reduction  in  the  overall  scatter. 

It  has  further  been  suggested  by  Sandham  and  Reynolds  (1989),  and  subse¬ 
quently  discussed  by  others,  e.g.,  Hall  et  al.  (1993),  that  Mc  does  not  provide 
a  universal  compressibility  measure,  even  within  a  set  of  experiments  performed 
in  a  given  facility.  Examples  of  this  non-universality  are  evident  in  Fig.  6.1b,  in 
reduced  values  of  8/8o  reported  by  Papamoschou  &  Roshko  (1988)  and  Hall  et 
al.  (1993).*  These  low-growth-rate  data  were  recorded  in  flows  characterized  by 
extreme  freestream  density  ratios,  e.g.,  s  ~  0.058  for  the  argon/helium  Case  9  flow 
of  Hall  et  al.  (1993).  Lu  &  Lele  (1994)  have  suggested  that  this  behavior  indicates 
that  the  departure  from  a  universal  collapse  is  a  systematic  deviation,  and  not  due 


* 


See,  by  way  of  example,  Case  9  of  Papamoschou  &  Roshko  (1988)  and  Cases  7-9  of  Hall  et 
al.  (1993) 
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Fig.  6.1a  Normalized  compressible  shear-layer  growth  rates  as  parameterized  by  the 
total  convective  Mach  number,  Mc  (Eq.  1.5).  Normalization  by  Eq.  1.2, 
with  CV s  as  suggested  (when  documented)  by  the  various  investigators. 

to  experimental  error.  This  systematic  departure  is  indicative  of  a  shortcoming  on 
one  of  two  counts: 

1.  An  inaccurate  representation  of  the  (normalizing)  incompressible  shear- 
layer  growth  rate,  Sq/x  (Eq.  1.2), 

or, 

2.  An  inaccurate  representation  of  flow  compressibility,  as  provided  by  the 
total  convective  Mach  number,  Mc  (Eq.  1.5). 


These  two  possibilities  will  be  explored  below. 
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FIG.  6.1b  Normalized  compressible  shear-layer  growth  rates  as  parameterized  by 
the  total  convective  Mach  number,  Mc  (Eq.  1.5).  Normalization  provided 
by  Eq.  1.2,  with  C$’ s  determined  as  described  in  text. 

6.2  The  incompressible-flow  growth  rate 

Lu  &  Lele  (1994)  recently  reported  an  improved  collapse  of  the  data  with 
the  use  of  a  modified  incompressible  growth-rate  estimate,  i.e.,  Option  1.  Their 
model  relied  on  an  incompressible-flow  linear-stability  analysis,  employing  the  most 
unstable  eigenvalue,  — (a,)0,  to  estimate  incompressible-flow  growth  rates. 

Compressible  shear-layer  growth  rates  appear  to  be  well-described  by  such  a 
correlation.  For  example,  Fig.  5  of  Lu  &  Lele  (1994)  illustrates  that, 

-a*  ~  -  ,  (6.1) 

x 

for  compressible-flow  conditions.  The  success  of  such  a  correlation  is  likely  due 
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to  the  relatively-small  spatial-growth  rates  and  reduced  turbulent-fluctuation  levels 
characteristic  of  compressible  turbulent  shear  layers,  features  which  enhance  the 
formal  validity  of  assumptions  inherent  to  linear-stability  analysis.  Interestingly, 
as  noted  by  Lu  &  Lele,  their  Fig.  5  also  suggests  a  rather  poor  correlation  for 
incompressible  flows  ( e.g .,  Cases  10  and  11  of  Hall  et  al.  1993). 


An  additional  finding  of  such  analyses,  e.g.,  Zhuang  (1990),  was  that  the  ratio 
of  the  compressible-flow  to  incompressible-flow  most-unstable  eigenvalue  was  found 
to  be  nearly  a  function  of  Mc,  for  a  wide  range  of  parameters,  i.e., 


-f—r-  -  fn(Mc)  , 
\&i  jo 


(6.2) 


which,  when  combined  with  Eq.  6.2  yields, 


8/x 
— (CKi)o 


~  fn(Mc)  . 


(6.3) 


The  near-functional  dependence  observed  by  Lu  Sz  Lele  (1994)  is  clear  from  Eq.  6.3. 
An  issue  remains,  however,  as  to  whether  the  incompressible-flow  growth  rate,  80/x, 
is  correlated  with  the  most-unstable  eigenvalue  from  an  incompressible-flow  linear- 
stability  analysis,  -(a,)0.  This  issue  is  especially  important  at  extreme  values  of 
the  freestream-density  ratio,  s ,  a  regime  in  which  linear-stability  analysis  has  not 
been  widely  tested,  as  will  be  discussed  below.  Unfortunately,  a  comparison  of 
their  model  with  existing  experimental  data,  e.g.,  Brown  Sz  Roshko  (1974),  was  not 
provided  by  Lu  &  Lele. 


Estimates  of  the  large-scale-structure  based  models  of  Brown  (1974)  and  Dimo- 
takis  (1986)  are  tested  with  data  from  Brown  &  Roshko  (1974)  in  Figs.  6.2.  These 
data  constitute  a  fairly  complete  survey  of  isolated  density-ratio  effects,  cf.  Eq.  1.3 
and  Ch.  3  above,  with  reported  values  for  80/x  at  density  ratios  of  s  =  1/7, 1,  &  7 
and  velocity  ratios  of  r  =  0, 1/7,  &  Over  this  range  of  conditions,  no  sys¬ 

tematic  deviations  of  the  model  estimates  from  the  data  are  discernible,  indicating 
that  Eq.  1.2,  and  its  temporal-growth  counterpart,  provide  an  adequate  estimate  of 
incompressible  shear-layer  growth  rates. 


**  Data  was  not  reported  for  the  r  =  0,  s  =  7  flow  condition. 
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Fig.  6.2a  Test  of  Brown  (1974)  model  predictions  for  (temporally-evolving)  incom¬ 
pressible  shear-layer  growth  with  the  data  of  Brown  &  Roshko  (1974). 

The  linear-stability  estimates  of  Lu  &  Lele  (1994)  are  provided  for  a  velocity 
ratio  of  r  =  1/2,  which  does  not  correspond  to  any  existing  set  of  experimental 
data.  As  a  result,  a  test  of  their  proposal  that  is  analogous  to  that  of  Figs.  6.2  is 
not  possible.  In  the  interest  of  evaluating  their  proposal,  however,  a  comparison  of 
the  three  estimates,  at  r  =  1/2,  is  presented  in  Fig.  6.3.  Not  surprisingly,  at  low 
density  ratios,  linear-stability-model  estimates  are  significantly  lower  than  those 
of  either  structure-based  model,  as  required  for  the  collapse  observed  by  Lu  & 
Lele.  This  reduction  is  extremely  sensitive  to  the  mean-profile  shape  imposed  upon 
the  linear  equations,  however,  as  a  hyperbolic-tangent  profile^  yields  estimates  at 

'  Lu  &;  Lele  employed  a  set  of  profiles  based  on  a  solution  to  the  laminar  thin  shear-layer 
equations,  with  the  two  freestreams  assumed  to  have  the  same  composition,  cf.  Hall  et  al.  (1993) 
Case  9  He/Ar  flow. 
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Fig.  6.2b  Test  of  Dimotakis  (1986)  model  predictions  for  (spatially-evolving)  incom¬ 
pressible  shear-layer  growth  with  the  data  of  Brown  &  Roshko  (1974). 

small  s  which  exceed  that  of  both  structure-based  models  (Zhuang  1995,  private 
communication).  In  addition,  a  non-monotonic  behavior  of  growth  rate  with  density 
ratio  is  indicated  for  the  linear-stability  estimates,  in  contrast  to  the  structure- 
based-model  estimates. 
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Fig.  6.3  Comparison  of  Brown  (1974),  Dimotakis  (1986),  and  Lu  &:  Lele  (1994) 
model  estimates  for  incompressible  shear-layer  growth  rates  at  r  =  1/2. 
Model  predictions  normalized  to  coincide  at  s  =  1. 


6.3  An  alternate  flow-compressibility  measure 

Turbulent  velocity  fluctuations  are  intimately  linked  to  shear-layer  growth.  Ex¬ 
perimental  support  for  such  an  association  includes  observation  of  decreased  (rms) 
velocity  fluctuations  with  reduced  shear-layer  growth  rates  by  Goebel  &  Dutton 
(1991),  in  addition  to  measurements  of  diminished  values  of  u'v'  at  compressible- 
flow  conditions  by  Samimy  &:  Elliott  (1990).  More  recently,  this  relationship  has 
been  formally  shown  in  terms  of  turbulent  kinetic  energy  by  Sarkar  (1995),  and 
in  terms  of  Reynolds-stress  production  by  Vreman  et  al.  (1996)  and  Freund  et 
al.  (1997). 


-  81- 


There  also  exists  considerable  evidence  to  suggest  that  growth-rate  suppression 
is  due  to  inviscid  energy-exchange  mechanisms,  and  not  due  to  increased  viscous 
(dilatational)  dissipation,  as  was  suggested  in  the  past  (e.g.,  Zeman  1990).  The 
data  indicate  no  dependence  of  growth  rate  on  Reynolds  number,  as  would  be 
expected  if  a  viscous  mechanism  were  dominant.  In  particular,  a  fortuitous  balance 
would  be  required  for  self-similar  (conical)  layer  growth,  as  the  Reynolds  number, 
Reg,  is  increasing  as  the  flow  evolves  downstream  (Dimotakis  1991a).  The  success 
of  inviscid  linear-stability  analyses  in  predicting  the  relative  reduction  in  growth 
rate,  e.g .,  Eq.  6.2,  lends  further  support  to  the  notion  that  the  dominant  growth- 
suppression  mechanisms  are  inviscid. 


In  the  most  general  terms,  flow  compressibility  acts  to  couple  kinetic  and  ther¬ 
mal  energy,  for  example,  in  the  steady-flow  kinetic-to-thermal  energy  (enthalpy) 


ratio, 


U1=hj-h=  J-1 

2  h  h  2 


(6.4) 


where  U  is  the  velocity,  and  M  the  Mach  number.  This  scaling  is  ubiquitous  in  all 
forms  of  the  dimensionless  energy  equation  (e.g.,  Lagerstrom  1964),  even  in  the  full 
unsteady  form.  Such  a  balance  must  be  observed  by  the  conversion  of  kinetic  energy 
into  thermal  energy  by  the  variety  of  mechanisms  available  to  compressible  turbu¬ 
lence.  Such  mechanisms,  e.g.,  dilatation,  are  unavailable  to  incompressible  flow 
and  must  observe  Eq.  6.4,  independently  of  the  specific  details  of  their  dynamics 
and  character.  Consequently,  as  flow  compressibility  increases,  velocity-field  fluctu¬ 
ations  can  no  longer  take  place  autonomously,  but  instead  result  in  corresponding 
thermal-energy  (enthalpy)  fluctuations. 


In  a  shear-layer,  the  freestream-density  ratio,  s,  is  related  to  the  speed-of-sound 


ratio,  i.e., 


(6.5) 


If  we  consider  a  shear  layer  with,  for  example,  a  small  density  ratio  (e.g.,  Case  9 
of  Hall  et  al.  1993),  a  correspondingly  large  speed-of-sound  ratio  is  indicated  by 
Eq.  6.5.  In  such  a  low-s  (or  conversely,  high-s)  flow,  however,  the  estimate  of  flow 
compressibility  given  by  Mc  (Eq.  1.5)  will  be  dominated  by  the  large  speed-of-sound 
value.  Given  the  propensity  of  these  flows  to  dynamically  attain  asymmetric  con- 
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vective  speeds,  resulting  in  correspondingly-large  convective-frame  Mach  numbers 
( e.g .,  Ch.  4),  such  an  estimate  is  likely  to  underestimate  the  flow  compressibility. 
Forming  the  group  suggested  by  Eq.  6.4,  using  the  (Galilean-invariant)  freestream 
velocity  difference,  AU  and  the  minimum  speed  of  sound,  yields  such  a  maximum 
compressiblity  scale, 

x  AU , 

or  explicitly  in  terms  of  the  freestream  values, 


(6.6a) 


1.2 


1.0 


0.8 


0.4 


0.2 


0.0 

0  12  3  4 

K 0 

Fig.  6.4  Normalized  (as  in  Fig.  6.1b)  compressible  shear-layer  growth  rates  param¬ 
eterized  by  the  proposed  compressibility  parameter  IIC  (Eqs.  6.6).  Solid 
curve  is  Eq.  6.7. 
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To  provide  an  estimate  of  the  relative  magnitude  of  heat  release  in  a  chemically¬ 
reacting  flow,  with  respect  to  these  measures,  it  is  useful  to  compare  the  change 
in  (bulk)  fluid  enthalpy  from  chemical  reaction,  A/irxn,  to  the  normalizing  quantity 
employed  in  Eqs.  6.6.  An  estimate  for  this  enthalpy  change  is  provided  by  the 
bulk-dilatation  correlation  of  Dimotakis  (1991a).  As  an  example,  for  Case  5.3b  flow 
conditions,  the  highest  heat-release  condition  in  this  work, 


\/ A hrxxi  X  max 


vTi 


a i 


~  0.1 


corroborating  the  validity  of  the  assumptions  of  low  heat  release  here.  At  very  high 
levels  of  heat  release,  however,  interesting  couplings  may  be  operative. 


The  experimentally-measured  growth  rates,  normalized  as  in  Fig.  6.1b,  are  plot¬ 
ted  against  IIC  in  Fig.  6.4.  These  data  observe  a  near-functional  dependence  on  this 
parameter,  to  a  better  extent  than  on  Mc.  The  solid  curve  through  the  data, 

~  (1  +  all2c)~\  a  ~  4,  /?~0.5,  (6.7) 

o o 

is  a  form  inspired  by  the  static-to-stagnation  enthalpy  ratio,  Eq.  6.4  above.  At 
high-compressibility  conditions,  this  functional  form  is  in  agreement  with  the  high- 
Mach-number  behavior,  i.e.,  1/M,  suggested  by  Brown  &  Roshko  (1974).  The  solid 
curve  in  Fig.  6.4  is  not  extended  to  high-compressibility  conditions,  as  additional 
energy-transfer  mechanisms  e.g.,  wave/bounding-wall  radiative  resonances  (Zhuang 
et  al.  1990),  may  well  play  a  role  in  this  regime. 

The  systematic  deviations  observed  in  flows  with  extreme  density /speed-of- 
sound  ratios  are  absent  when  scaled  in  this  fashion.  Considering  that  many  flows  re¬ 
ported  in  the  literature  were  generated  with  thermodynamically-similar  freestream 
fluids,  e.g.,  the  air/air  experiments  of  Chinzei  et  al.  (1986),  resulting  in  near-unity 
values  of  02  /a\  and  uniform  values  of  7,  it  is  not  surprising  that  the  previous  col¬ 
lapse  of  normalized  growth-rate  data  was  viewed  as  adequate  by  many  investigators. 
For  such  conditions, 

a2  =  01  and  72  =  7i  =  7  =£■  nc  =  2^/7  —  1  x  Mc  ,  (6.8) 

we  would  expect  the  limited  (Mc)  collapse  to  be  preserved  by  the  parameterization 
in  terms  of  IIC,  as  is  borne  out  by  Fig.  6.4. 
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6.4  Conclusions 


An  alternate  shear-layer  flow  compressibility  parameter  was  motivated  on  the 
basis  of  shear-layer  growth-rate  data  and  a  view  of  compressibility  as  an  energy- 
transformation  mechanism.  This  parameter, 


nc 


max 

2 


-  1 


a; 


x  A  U, 


which  differs  from  the  commonly-utilized  total  convective  Mach  number,  Mc,  at  ex¬ 
treme  freestream-density  ratios,  provides  a  good  collapse  of  experimentally-observed 
growth  rates.  In  addition,  a  functional  form  based  on  the  static-to-stagnation  en¬ 
thalpy  ratio, 

7-  —  (l  +  all*)  13  ,  a  ~  4,  j5  ~  0.5  , 
o0 

provides  a  good  representation  of  the  experimental  data,  for  IIC  <  3. 


With  regard  to  some  envisioned  supersonic-combustion  scenarios,  e.g.,  H2-a.ir 
propulsion,  in  which  extreme  density  ratios  and  associated  speed-of-sound  ratios 
will  doubtlessly  be  encountered,  this  parameterization  may  prove  useful. 
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CHAPTER  7 

Conclusions 


Experimental  investigations  of  several  aspects  of  high-Reynolds-number,  tur¬ 
bulent,  shear-layer  flows  were  undertaken  in  the  GALCIT  Supersonic  Shear  Layer 
Facility.  These  consisted  of  both  non-reacting  and  chemically-reacting  flows  and 
yielded  quantitative  measurements  of  large-scale  dynamics  as  well  as  small-scale 
(molecular)  mixing,  with  chemically-reacting  flows  employing  the  (H2  +  N0)/F2 
chemical  system  in  the  fast-kinetic  regime.  The  properties  of  this  canonical  tur¬ 
bulent  flow  are  known  to  depend  on  a  host  of  parameters.  This  work  focused  on 
effects  of  inflow/initial  conditions,  effects  of  freestream  boundary  conditions  (super¬ 
sonic/subsonic  flow),  and  effects  of  compressibility. 

Individual  conclusions  summarizing  the  findings  of  each  part  of  the  study  are 
found  following  each  Chapter.  A  summary  is  offered  below. 

Inflow  conditions  —  Effects  of  perturbed  inflow  conditions,  i.e.,  tripping  the 
boundary  layer  flows  on  the  splitter  plate  upstream  of  shear-layer  formation,  were 
observed  to  persist  for  distances  that  must  be  considered  far  downstream,  in  a  high- 
Re  ( Res  ~2x  105),  incompressible  flow.  Growth  rates  and  mixed-fluid  properties 
were  significantly  altered  by  tripping  the  high-speed  splitter-plate  boundary  layer. 
The  mixed-fluid  (scalar)  field,  at  a  downstream  distance  x/9\  ~  3300  and  a  “pairing 
parameter”  of  47,  changed  from  a  non-marching  probability- density  function  (pdf), 
to  a  marching  pdf,  when  the  boundary  layer  was  tripped.  Several  other  perturbation 
configurations  were  explored,  e.g.,  a  low-speed-side  trip,  with  little  observed  effect 
on  either  growth  or  mixing.  These  observations  are  indicative  of  shear-layer  behavior 
that  depends  on  not  only  local-flow  properties,  but  also  on  upstream  conditions, 
reminiscent  of  results  from  low-dimensionality,  non-linear  (chaotic)  systems. 
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Freestream  conditions  —  The  presence,  or  absence,  of  subsonic  (laboratory- 
frame)  regions  on  the  flow  was  explored  in  a  set  of  experiments  where  the  freestreams 
were  bi-subsonic,  bi-supersonic,  and  (hybrid-flow)  supersonic/subsonic.  Fast-mode 
turbulent-structure  convection  speeds  were  observed  in  both  bi-subsonic  and  bi- 
supersonic  flows,  at  variance  with  an  existing  empirical  rule  for  mode  selection.  This 
rule  had  suggested  that  the  presence  of  (one)  subsonic  freestream  was  responsible 
for  fast-mode  dominance,  i.  e.,  Mc\  <C  Mc 2.  Growth  rates  were  suppressed  with 
increasing  flow  compressibility  (Mc),  in  accord  with  previous  data.  This  further 
indicates  that  turbulent  shear  layers  are  relatively  insensitive,  at  least  for  the  flows 
investigated,  to  laboratory-frame  subsonic  regions.  In  the  transonic-flow  regime, 
however,  global  flowfields,  i.e.,  the  two  freestreams  coupled  by  the  shear  layer,  were 
found  to  be  sensitive  to  imposed  boundary  conditions. 

Compressibility  and  mixing  —  Effects  of  compressibility  on  molecular  mix¬ 
ing  were  investigated  with  chemically-reacting  flip  experiments,  at  two  different 
freestream  flow  conditions.  These  conditions  were  characterized  by  nominally-fixed 
freestream-velocity  and  -density  ratios,  local-flow  Reynolds  numbers,  and  reac¬ 
tant  compositions,  but,  a  changing  convective  Mach  number,  from  Mc  ~  0.25  to 
Mc  ~  0.47.  This  feature  allows  differences  in  the  measurements  to  be  attributed 
to  the  change  in  compressibility  alone.  Large-scale  growth  rates  were  reduced  with 
increasing  compressibility,  in  accord  with  previous  findings.  Mixed-fluid  properties 
exhibit  several  changes  associated  with  increasing  compressibility.  In  particular, 
compressibility  has  little  effect  on  the  H2-rich  part  of  the  “flip”,  but  increases  the 
peak  temperature  rise  of  the  F2-rich  part  by  approximately  25%,  yielding  a  com¬ 
mensurate  change  in  mixed-fluid  composition  ratio.  Only  a  small  change  in  the 
mixed-fluid  fraction  was  found  to  be  associated  with  the  increase  in  compressibility. 

Compressibility  and  dynamics  —  An  alternative  shear-layer  flow  com¬ 
pressibility  parameter,  derived  from  consideration  of  compressibility  as  an  energy- 
transformation  mechanism,  was  proposed.  A  test  of  incompressible  growth-rate 
models  was  conducted,  indicating  that  large-scale-structure  models  provide  an  ac¬ 
curate  representation  of  incompressible-flow  data  and  are  adequate  for  normaliza¬ 
tion  of  compressible-flow  data.  Normalized  shear-layer  growth-rate  data  are  well 
represented  by  the  proposed  parameter. 
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APPENDIX  A 


Color-schlieren  imaging 


A.l  Introduction 

The  turbulent  shear  layer  provides  an  ideal  subject  for  flow  visualization  by 
refractive-index  sampling  techniques,  e.g .,  schlieren  photography,  because  of  the 
substantial  gradients  generated  by  vigorous  mixing  processes,  as  well  as  those  pro¬ 
duced  by  dilatation  in  compressible  flows.  Such  visualizations,  employing  conven¬ 
tional  black-and-white  schlieren  and  shadowgraph  techniques,  have  provided  sig¬ 
nificant  information  regarding  structure  and  large-scale  dynamics  in  these  flows, 
from  the  identification  of  spanwise- coherent  large-scale  structures  (Brown  &  Roshko 
1974),  to  the  apparent  absence  of  such  features  in  compressible  flows  (e.g.,  Goebel 
&  Dutton  1991,  Hall  et  al.  1993,  and  Clemens  &  Mungal  1995). 

The  information  content  of  such  black-and-white  visualizations,  however,  can 
only  record  a  scalar  contraction  of  the  two-dimensional  (spanwise-integrated)  index- 
of-refraction  gradient  field.  A  black-and-white  schlieren  system  typically  represents 
a  single  component  of  the  field,  corresponding  to  the  direction  perpendicular  to 
the  knife  edge.  These  one-dimensional  data  can  often  lead  to  difficulties  in  inter¬ 
preting  flow  structure,  as  seen  in  the  black-and-white  schlieren  data  (Figs.  6.7  and 
6.8)  of  Hall  (1991),  in  which  a  combination  of  (non-simultaneous)  vertical  and  hor¬ 
izontal  knife-edge  realizations  was  required  to  observe  the  details  of  the  unsteady 
pressure- wave  system  in  a  confined,  supersonic,  shear-layer  flow.  As  discussed  by 
many  authors,  most  recently  in  reviews  by  Settles  (1985)  and  Kleine  &  Gronig 
(1991),  the  additional  information  in  a  color  realization  allows  simultaneous  resolu¬ 
tion  of  both  components  of  the  spanwise-integrated,  refractive-index-gradient  field. 
In  interpreting  complex,  multi- dimensional  fields  characteristic  of  high-speed  mixing 
environments,  the  additional  data  provided  by  color  can  be  of  significant  benefit. 
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Color  schlieren  techniques  are  broadly  classified  into  two  categories:  the  color- 
dissection  technique  of  Cords  (1968)  and  the  filtered-cutoff  technique  of  Rheinberg 
(1896).  The  former  system  employs  a  filter  mask  at  the  optical  source  plane  to 
“dissect”  the  light  into  different  color  components,  which  are  then  “recombined” 
at  the  observation  plane.  Application  of  this  technique  has  found  (relatively)  re¬ 
cent  use  in  high-speed  flows,  e.g.,  Kleine  &  Gronig  (1991).  Although  it  alleviates 
the  diffraction-induced  resolution  degradation  of  a  slit-like  cutoff-filter  arrangement, 
e.g .,  Rheinberg  (1896),  a  serious  drawback  of  the  dissection  technique  is  the  rela¬ 
tively  severe  illumination  requirement,  since  a  large  fraction  of  the  light  incident  on 
the  source  mask  is  blocked.  The  configuration  discussed  by  Kleine  &  Gronig,  for 
example,  required  a  custom-built  spark  source  significantly  more  energetic  than  was 
required  by  their  equivalent  (same  optical  components)  black-and-white  system.  ^ 

The  latter  system,  a  variant  of  which  is  employed  here,  is  a  direct  extension  of 
classical  monochrome  schlieren  techniques,  and  as  such,  potentially  shares  similar 
sensitivity,  dynamic  range,  and  illumination  demands.  The  only  modification  to 
the  black-and-white  system  is  the  introduction  of  a  color-filter  array,  in  place  of 
the  usual  knife  edge,  at  the  stop  plane.  It  is  often  stated,  e.g.,  Settles  (1985),  that 
such  a  system  suffers  from  a  degradation  in  spatial  resolution,  because  of  diffrac¬ 
tion  at  the  stop  plane.  In  the  case  of  “colored-slit”  or  “lattice-filter”  stop-plane 
arrays,  diffraction  resulting  from  the  arrangement  of  (effective)  slits  can  certainly 
be  significant,  especially  in  a  system  of  reasonable  physical  dimensions.  If,  however, 
an  array  which  codes  color  in  the  same  manner  as  a  knife  edge  codes  intensity  is 
employed  (see  Fig.  A.l),  then  the  deleterious  diffractive  effects  are  eliminated.  In 
such  a  system,  the  introduction  of  the  additional  dimension  provided  by  color  is  ac¬ 
companied  by  the  spatial  resolution  of  the  black-and-white  system,  at  a  comparable 
illumination  level.  This  will  be  demonstrated  below. 


t 


This  increase  in  illumination  energy  was  partially  achieved  by  extending  the  spark  duration, 
resulting  in  a  reduction  in  the  temporal  resolution. 
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A.2  Visualization  of  a  high-speed  turbulent  shear-layer  flow 

A  schematic  of  the  optical  system  employed  here  is  found  in  Fig.  A.l.  This 
“folded”  layout  is  based  on  a  design  due  to  Liepmann  (1947)  (see  discussion  by 
Herman  son  1985  and  Hall  1991),  which,  because  of  the  symmetrical  layout,  is  free 
of  coma.  The  astigmatism  in  the  system  is  minimized  via  the  use  of  large  focal- 
length  schlieren  heads  (/  ~  2m  first-surface  mirrors).  Black-and-white  visualiza¬ 
tions  of  shear-layer  flows,  recorded  with  this  apparatus,  possess  adequate  contrast 
and  dynamic  range,  e.g.,  Hall  et  al.  (1993),  and  as  a  consequence,  extension  to  a 
filtered-cutoff  system  would  be  expected  to  yield  good  results.  This  performance  is 
primarily  the  result  of  an  optically-long  integration  path  (6  =  15.3  cm  test-section 
span,  maintained  at  atmospheric  static  pressure)  which  enhances  the  magnitude  of 
recorded  gradients. 

The  color  cutoff  mask  was  fabricated  from  120°-wedge  segments  of  red,  green, 
and  blue  (Lee  #’s  019,  124,  and  183,  respectively)  filters. 11  As  indicated  above,  this 
extension  to  a  cutoff-mask  color  configuration,  in  contrast  to  a  dissection  technique, 
does  not  require  a  modification  of  the  illumination  source.  In  particular,  we  axe  still 
able  to  employ  a  commercially-available  Nanospark  N-789B  spark  source  driven  by 
a  N-487B  high-voltage  source  (20  ns  FWHM  illumination  duration),  fast  enough  to 
temporally  resolve  high-speed  flows.  If  an  equivalent-sensitivity  dissection  technique 
had  been  employed,  the  development  of  a  longer-duration,  higher-output  source 
would  have  been  required,  with  a  probable  loss  of  temporal  resolution. 

An  image  of  a  /  =  10  m  lens,  as  viewed  through  the  optical  system,  is  re¬ 
produced  in  the  schematic  (Fig.  A.l).  This  lens  produces  a  known  (light-beam) 
deflection,  analogous  to  deflections  produced  by  flow-induced  index-of-refraction 
gradients,  providing  a  calibration  reference.  The  maximum  light-ray  deflection, 
observed  at  the  edge  of  the  lens,  is  estimated  to  be  approximately  1  milliradian. 

A  composite  image  of  a  subsonic,  but  moderate-to-high  speed,  shear-layer 
flow  (Case  5.2a  conditions),  recorded  with  the  system  described  above,  is  shown 


#  Other  masks  were  fabricated  from  Kodak  Wratten  #’s  25  (red),  58  (green),  47  (blue),  and  8 
(yellow). 
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Fig.  A.l  Schematic  of  schlieren  photography  system.  Circular  pattern  at  upper  left 
is  an  /  =  10  m  lens,  as  imaged  through  the  schlieren  optics,  to  illustrate 
the  performance  of  the  system. 


Fig.  A. 2  Composite  color  schlieren  photograph  Case  5.2a  flow. 


in  Fig.  A. 2.  As  with  all  of  the  images  presented  in  this  work,  the  image  consists  of 
two  uncorrelated  photographs,  taken  during  different  runs,  to  visualize  both  the  up- 
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stream  and  downstream  regions  of  the  flow.  In  this  flow,  the  freestream  gas  mixtures 
have  been  tailored  so  that  freestream  densities  are  nearly  uniform  (s  ~  1),  with  the 
sckliere  largely  produced  by  the  temperature  rise  from  the  chemical  reaction. 

In  the  flow  image,  the  large-scale,  spanwise-coherent,  vortical  structures,  typical 
of  these  flows,  are  clearly  visible.  Compared  to  black-and-white  visualizations  of 
such  flows,  e.g.,  Brown  &  Roshko  (1974),  we  can  glean  additional  information  from 
this  realization,  because  of  the  additional  dimension  of  index-of-refraction  gradient 
data  recorded  in  the  color  visualization.  For  example,  the  local  (instantaneous) 
structure  of  the  temperature-rise  field,  in  this  chemically-reacting  flow,  is  seen  to 
be  nearly  symmetric  with  respect  to  the  large-scale  structures,  an  observation  that 
cannot  be  made  on  the  basis  of  monochrome  schlieren  data  without  assumptions  as 
to  local-gradient  directions.  Also  visible  in  this  image  is  the  top-guidewall  boundary 
layer,  which  supports  a  temperature  gradient  between  the  cooler,  high-speed  flow 
(Tflow  ~  265  K)  and  the  warmer,  room-temperature  (Twaii  —  295  K)  guidewall. 

Several  general  inferences,  with  regard  to  performance  of  this  color- schlieren 
system,  are  also  obtained  from  this  image.  First  and  foremost,  it  is  clear  from  the 
small-scale  detail  present,  e.g.,  the  upper-guidewall  boundary-layer  structure,  that 
the  spatial  resolution  is  relatively  high  and  has  not  been  significantly  degraded  by 
diffractive  effects.  Of  course,  this  is  not  surprising,  since  with  regard  to  diffractive 
effects,  the  filter  arrangement  employed  is  equivalent  to  the  classical  knife-edge. 
Second,  the  high  spatial  resolution  has  been  achieved  with  efficient  use  of  the  illu¬ 
mination  source,  relative  to  a  dissection-based  system.  Consequently,  the  high  tem¬ 
poral  resolution  associated  with  short-duration  (20  ns)  spark-source  illumination  is 
retained.  It  is  also  evident  from  the  range  of  colors  generated  by  the  calibration  lens 
that  the  sensitivity  and  dynamic  range  of  this  system  are  qualitatively  comparable, 
if  slightly  reduced,  to  that  of  a  similar  dissection-based  system  (c/.,  Kleine  &  Gronig 
1991),  at  a  much-reduced  illumination  level. 


-  92- 


APPENDIX  B 

Compressible  shear-layer  growth-rate  normalization 


As  discussed  in  Ch.  6,  values  employed  for  Cg  (Eq.  1.2)  in  the  normalization 
of  compressible  shear-layer  growth-rate  data  are  usually  not  measured,  but  in¬ 
stead  chosen  to  match  incompressible-flow  data,  t.g.,  of  Brown  &;  Roshko  (1974). 
Given  the  variance  exhibited  in  Eq.  1.3,  a  range  which  is  likely  also  exhibited  in 
compressible-flow  data,  an  attempt  was  made  to  reconcile  the  entire  data  set  from 
various  facilities/experiments  by  recomputing  the  C^’s  as  follows. 

Since  most  individual  data  sets  contain  measurements  near  Mc  —  0.5,  the 
empirical  correlation  suggested  by  Dimotakis  (1991a), 

4-  ~  /(Mc)  =  0.8e-3M'  +o.2  ,  (B.l) 

c>o 

was  employed  to  as  a  local  interpolator  of  the  normalized  growth  rate,  6 /So,  for  a 
particular  data  point  at  Mc,re f,  close  to  Mc  =  0.5.  This  correlation  was  based  on 
the  data  available  at  the  time,  and  uses  a  functional  form  similar  to  that  employed 
by  Zhuang  (1990)  in  an  investigation  of  compressible-flow  linear-stability  analysis. 
It  has  a  different  functional  form  (exponential  vs.  power  law)  than  the  correlation 
offered  in  Ch.  6,  which  uses  a  different  compressibility  parameterization. 

This  estimate  was  then  used  to  calculate  Cg,  i.e., 

c>-  Qr6it  >  (B-2) 

where  8o/(xCg)  is  computed  using  Eq.  1.2,  ^(MCjre f)  is  the  growth  rate  measured 
at  Mc  =  MC)ref,  and  /(Mc,ref)  is  computed  using  Eq.  B.l.  Several  comments  are 
warranted: 

1.  The  data  of  Sirieix  &:  Solignac  (1966)  are  normalized  at  MC)ref  =  1.125, 
the  lowest-compressibility  flow  reported. 
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2.  The  (single-point)  data  of  Ikawa  &  Kubota  (1975)  employs  the  same  Cs 
as  that  of  Samimy  &  Elliott  (1990). 

3.  No  correction  is  made  for  the  (single-point)  data  of  Barre  et  al.  (1994). 

This  procedure  results  in  the  following  values  for  Cs,  where  Cs,0  are  values 
suggested  (or  deduced  from)  in  the  original  publications. 


Sirieix  &  Solignac1: 
Ikawa  &  Kubota2: 
Chinzei  et  al.: 
Papamoschou  &:  Roshko: 
Samimy  &  Elliott: 
Goebel  &  Dutton: 
Hall  et  al.: 
Barre  et  al.3: 
Clemens  &  Mungal: 

Table  C.l:  Details  of  Cs 


-17c>ref 

c,ref  ) 

1  So 
Cg  x 

1.13 

0.218 

1.50 

0.991 

0.242 

1.59 

0.591 

0.481 

0.805 

0.555 

0.517 

0.400 

0.517 

0.559 

0.881 

0.453 

0.632 

0.570 

0.511 

0.565 

1.16 

0.623 

0.450 

1.19 

0.493 

0.586 

1.05 

determination.  Superscripts 


Cs 

Cfi)0 

0.0345 

0.106 

0.140 

0.064 

0.189 

0.105 

0.045 

0.116 

0.140 

0.026 

0.132 

0.140 

0.093 

0.189 

0.105 

0.038 

0.105 

0.0825 

0.114 

0.174 

0.170 

0.0338 

0.0606 

0.0606 

0.100 

0.163 

0.17 

refer  to  entries  in  list  above. 


The  present  work  employs  the  same  values  as  Hall  et  al.  (1993),  as  these  exper¬ 
iments  were  in  the  same  facility,  and  employ  the  same  (visual)  growth-rate  measure. 
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